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POSSIBILITY OF EXCHANGE MAGNETOSTRICTION 
YIELDING NEGATIVE THERMAL EXPANSION IN SOLID He® 


D. S. Rodbell 
General Electric Research Laboratory, Schenectady, New York 
(Received April 10, 1961; revised manuscript received May 24, 1961) 


According to Primakoff! the nuclear spins of 
solid He® should be antiferromagnetically cou- 
pled with a transition temperature that decreases 
with increasing pressure. The essential features 
of antiferromagnetism have been observed by 
Meyer, Adams, and Fairbank? to occur in the 
region of the He® phase diagram near 0.1°K for 
P>70 atm. This observation is within order-of- 
magnitude agreement with the predictions of Ber- 
nardes and Primakoff** and lends support to the 
exchange interactions so far calculated.’»** This 
exchange is believed to be a steep function of lat- 
tice parameter®® in the highly compressible (K 
~5x107® cm?/dyne) He® lattice. The purpose of 
this note is to point out that the existence of (a) an 
exchange interaction dependent upon interatomic 
separation and (b) a highly compressible lattice 
are sufficient conditions to give rise to negative 
thermal expansion behavior. Such behavior has 
been observed® in certain regions of the solid He® 
phase diagram. We outline the features of this 
mechanism’ below. Bernardes and Primakoff® 
have shown such behavior to be contained in their 
spin-wave treatment of the He*® problem; we believe 
that the phenomenological approach used here pro- 
vides some additional physical insight. The phys- 
ical basis of the mechanism is provided by the 
fact that the larger the exchange interaction, the 
lower the system’s free energy. Thus a sponta- 
neous distortion in the direction of increasing the 
exchange will lower the energy but must be paid 
for with the strain energy of the distorted lattice. 
In addition a higher effective Curie temperature 
provides more magnetization at a given tempera- 


ture and hence lowers the spin entropy. 
We now write the Gibbs function (G =U -TS+PV) 
for a spin $ system and for a volume v, as 





-v.\?71 
= -H = 2 1 v Vo be 
G. M .ocos($/2) NRT 0 cos¢ + i( - ) 


0 o /K 
field exchange strain 





v-v 
-T 2). 
Sn +5 stew +P( Vo ) (1) 
entropy pressure 


In this expression, H is the magnetic field, M, 
the sublattice saturation magnetization, and o the 
relative magnetization at temperature T; ¢ is the 
angle between the assumed sublattices, being 0 
for ferromagnetism and 7 for antiferromagnetism. 
N is the number of magnetic atoms per unit vol- 
ume, k is Boltzmann’s constant, v is the specific 
volume, v, is the specific volume in the absence 
of a magnetic interaction, and K is the compres- 
sibility. The Curie temperature T.. which we use 
to measure the exchange interaction in the molec- 
ular field approximation is assumed to be given 
by 


T¢ =T,{1+B[( -v,)/vo}}, (2) 


where T, is the “undistorted” Curie temperature 
(that may be positive or negative corresponding 
to ferromagnetism or antiferromagnetism) and 
6 measures the slope of the dependence of T/T, 
with specific volume (and may be either positive 
or negative). 

The free energy depends upon the lattice thermal 
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expansion coefficient, a p as this parameter meas- 
ures the change of lattice entropy with volume. We 
assume that Sj, +4 =S,+(a@,/K)(v -v,), an expression 
that is justifiable for a constant a, but otherwise 
is simply assumed for this treatment. The spin 
entropy does not enter directly into the present 
determination of expression (3) below; it is im- 
portant in determining o(T), however, and there- 
fore determines the temperature dependence of (3). 
We have found situations® where the spin entropy 
may dominate a phase transition. The fact that 
the spin system entropy may play a dominant role 
has been emphasized by Goldstein® and Bernardes 
and Primakoff** with regard to solid He’. 

The specific volume that minimizes the free 
energy (1) is found to be 


U-Uo 


(3) 





=a T - PK + 4NRKT,Bo’ cos. 


Vo : 
equil 


We take 7, negative, 8B positive, and in the ab- 
sence of an external field, ¢=7 as descriptive 

of the antiferromagnetic state of He*.*~> The de- 
pendence of (3) on temperature is the thermal 
expansion to be expected in such a system, and 

in Fig. 1 we indicate the qualitative features of 

(3) where to the “normal” thermal expansion must 
be added a term proportional to o”. This term 
disappears with the magnetization as temperature 
increases and results in the possibility of a change 
in sign of the thermal expansion coefficient. The 
necessary condition for negative thermal expan- 
sion behavior is that a,< §NKkf be true.*® This 

is indeed the case for solid He* using the measured 


TOTAL THERMAL EXPANSION 
WITH SLOPE CHANGE 


SPECIFIC 
VOLUME 





Pg 
NOR MAL THERMAL EXPANSION 


values of (a,,K)**™ and with (1/T ,)(8T ,/@P) ~ 107 ’ 
from Bernardes and Primakoff® which corresponds | 
to B =2. 

If the exchange energy vs lattice parameter has 
a monotonic dependence, then the distortion that 
in the antiferromagnetic state favors an expanded 
lattice will in the ferromagnetic state become a 
contraction, since this provides the largest abso- 
lute values of T,, (hence the lowest free energy) 
in each case. These views would argue in favor 
of a denser ferromagnetic solid than an antiferro- | 
magnetic one and would suggest that the exchange [ 
energy and its dependence upon atomic separation 
directly contributes to the phase change that yields) 
a denser solid at high pressure, as has been ob- 
served. The latter conclusion supports recent 
suggestions by Bernardes and Primakoff.* The 
detailed treatment® of the problem outlined here 
also evaluates the “magnetic phase diagrams” of | 
such systems and finds possible first- and second- 
order phase transitions occurring as pressure or 
temperature are varied. The predictions by Ber- 
nardes’ and Primakoff® of a transition from anti- 
ferromagnetism to ferromagnetism’ with an ex- 
ternal pressure of 200 atmospheres or greater 
(near 0.1°K) is in accord with our own “model” 
using Bernardes’ and Primakoff’s value for §; 
in addition we would also predict a range of tem- 
peratures where transitions from antiferromag- 
netism ~ paramagnetism ~ ferromagnetism would 
occur with increasing pressure. This last feature 
may have bearing upon the unusual susceptibility 
behavior that has been observed.” 

In conclusion, it should be noted that a depend- 






FIG. 1. The total thermal 
expansion of a highly com- 
pressible lattice with exchange 
interaction as described in the 
text. 
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ence of exchange on lattice parameter of opposite 
sign (8 negative) would not yield (for an antiferro- 
magnetic ground state) negative thermal expansion 


behavior since the exchange magnetostriction term 


would give rise to a contraction at T=0. It may 


thus be that the observation of negative thermal ex- 


pansion in solid He® is auxiliary evidence for the 
supposed dependence of exchange upon lattice pa- 
rameter. 

I have benefited from communications with 
Professor H. Primakoff on the subject of this 
note, the kernel of which originated in a more 
general problem on magnetic phase transitions 
done in conjunction with C. P. Bean. 





'H. Primakoff, Bull. Am. Phys. Soc. 2, 63 (1957). 

2H. Meyer, E. D. Adams, and W. M. Fairbank, 
Proceedings of the Seventh International Conference on 
Low-Temperature Physics, September, 1960 (Univer- 
sity of Toronto Press, Toronto, 1960), p. 339. 

3N. Bernardes and H. Primakoff, Phys. Rev. Letters 
2, 290 (1958). 

‘N. Bernardes and H. Primakoff, Phys. Rev. Letters 
3, 144 (1959). 

5N. Bernardes and H. Primakoff, Phys. Rev. 119, 
968 (1960). 

SE, Grilly, S. Sydoriak, and R. Mills, in Helium 3, 
edited by J. G. Daunt (Ohio State University Press, 
Columbus, Ohio, 1960). 








‘Cc, P. Bean and D. S. Rodbell, Bull. Am. Phys. 
Soc. 6, 159 (1961). 
8c, P. Bean and D. S. Rodbell (to be published). 

°L. Goldstein, Ann. Phys. 8, 390 (1959). 

This condition results from evaluating |odo/dT| at 
its maximum €~ Ofrom the Brillouin function for j=1/2 
where limg-»9|od0/aT |=3/2T). The effect of modifying 
the Brillouin function as occurs in this treatment will 
raise this value; we have therefore evaluated the lower 
bound of this exchange magnetostriction. 

''Estimating a7 from the high-temperature measure- 
ments to be less than 10~*/°*K while K ~5 x107° 
cm?/dyne. 

12n, Bernardes, in Helium 3, edited by J. G. Daunt 
(Ohio State University Press, Columbus, Ohio, 1960). 

134 transition of this type may be appreciably af- 
fected by externally applied magnetic fields and may 
be understood through the magnetic analog of the 
Clausius-Clapeyron relation, e.g., AT/T=HAM/L; 
for latent heat, L, a transition temperature change 
AT will result for an applied field H; AM is the dif- 
ference in net magnetization per unit volume between 
the ferromagnetic and antiferromagnetic states. If 
the latent heat is largely due to changes in the spin 
system entropy, then L may be small and AT large. 

“Note added in proof. M. J. Klein and R. D. Moun- 
tain [Phys. Rev. Letters 5, 363 (1960)] have pointed 
out that the thermal expansion must change sign when 
the transition from antiferromagnetism to ferromag- 
netism occurs. This is also evident in our Eq. (3) 
when ¢ changes from 7 to 0. 








FLUX AND ENERGY SPECTRA OF THE PROTONS IN THE INNER VAN ALLEN BELT 


J. E. Naugle and D. A. Kniffen 


National Aeronautics and Space Administration, Goddard Space Flight Center, Greenbelt, Maryland 
(Received June 8, 1961) 


In this experiment, for the first time, the flux 
and energy spectrum of the trapped protons in 
the inner Van Allen belt were measured as a 
function of position in the belt. Previous meas- 
urements!” of the proton energy spectrum have 
been made in emulsions which integrated over 
the entire flight paths of intercontinental ballis- 
tic missiles flown from Cape Canaveral. 

On September 19, 1960, a research rocket 
with a special nose cone developed for the ex- 
posure and recovery of emulsions was flown 
along a trajectory which closely corresponded 
to a magnetic meridian. The rocket payload 
was designed to permit a four-inch cylinder 
to be extended out through the front end of the 
nose cone, thereby reducing the amount of ma- 
terial surrounding the emulsion container to a 


minimum. Time resolution of the spectrum was 
obtained by exposing sections of a cylindrical 
emulsion stack to the ambient radiation through 
a port in the tungsten emulsion container. Fig- 
ure 1 shows the experimental arrangement. A 
point on the periphery of the stack was behind 
the port for approximately 80 seconds. During 
the rest of the flight that point was shielded by 
30.6 g/cm? of tungsten. Therefore, the track 
population at such a point is made up of particles 
which came in through the port during the 80- 
second exposure and a background of particles 
which passed through the tungsten during the 
remaining 900 seconds the payload was in flight. 
A portion of the stack which did not pass behind 
the port has been used to measure the background. 
Protons of energy 2 8 Mev were able to penetrate 
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FIG, 1. 25 circular G-5 


Film position at apogee. 
‘and G-0 nuclear emulsions were contained in a tung- 


sten cassette. The numbered points around the per- 
iphery of the stack refer to the sections of the emul- 
sion which are being analyzed. 


the aluminum port cover and be detected, where- 
as only protons with E 2 145 Mev were recorded 
in the emulsion after passing through the tung- 
sten shield. 

Figure 2 shows the trajectory of the rocket 
payload. The numbered points along the trajec- 
tory indicate the locations at which the corre- 
sponding points on the stack were behind the 
port. Table I gives the pertinent data regarding 
the three points 3, 4, 5, which have been ana- 
lyzed to date. 

A “line scan” was made around the periphery 
of the emulsions at points 3, 4, and 5. The scan 








was made approximately 500 microns from the 
edge of the emulsions. It was not possible to 
scan closer to the edge of the emulsion because 
of blackening. All particles which, at the scan 
line, were within +25° of a radius vector and 

+45° of the plane of the emulsion were recorded. 

The portion of the emulsion which was not ex- 
posed behind the port was also scanned, with the 
same criteria, to obtain an integral flux averaged 
over the entire flight and also to obtain a back- 
ground correction for the other sections of the 
stack. 

In order to obtain the differential energy spec- 
trum, all particles in the scans which had a pos- 
sible residual range in the emulsion, at the scan 
line, of at least 2.5 mm were followed until they 
either ended, interacted, or left the stack. The 
energies and mass of these particles were then 
determined from range, ionization, and scatter- 
ing measurements. 

In the analysis to date, no evidence has been 
found for particles other than protons crossing 
the scan line. However, the mass measurements 
do not rule out a small component of singly- 
charged particles heavier than protons. The 
thickness of the blackened edge of the emulsion 
varies around the periphery of the stack and is 
consistent with electrons of energy up to about 
600 kev stopping in the emulsion. There isa 
high background of grains for another 500 mi- 
crons in from the black edge, due possibly to 
the presence of soft x rays. 

The number of protons with energy greater 
than 100 Mev crossing the scan lines was ap- 
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Table I. Summary of results. 





Point 3: 


1600 km J(E=31) = 70+ 15 protons/cm’ sr sec 
28.4°N lat. 4n J(E 231) =900 +200 protons/cm? sec 
237.7°'W long. JE >31)* = 220 protons/cm’ sec 
B=0.231 gauss 
I=1.67 - 


J(E) =(3.2+0.6) x 108 /E*5 **5 protons/cm’ sec sr Mev; 10<E <50 Mev 


Point 4: 
1884 km J(E =31) = 120+ 20 protons/cm’ sr sec 
25.2°N lat. 4nJ(E =31)_ =1500 + 250 protons/cm? sec 
236 .5°W long. JE 231) =1000 protons/cm? sec 
B=0.198 gauss 
I=1.34r7 
e 
J(E) =(6.8 +0.3) x 10° /E*5 **? protons/cm’ sec sr Mev; 10<E <50 Mev 
Point 5: 
1600 km J(E =31) = 190+ 40 protons/cm’ sr sec 
21.5°N lat. 4nJ(E =31 Mev), = 2400 +500 protons/cm? sec 
235.4°W long. Jo(E =31 Mev) =1290 protons/cm? sec 
B=0.209 gauss 
1=0.946 .. 


J(E) =(1 +0.5) <10°/E*"*%3 protons/cm? sec sr Mev; 40<E < 100 Mev 
Integral flux behind 30.6 g/cm? of tungsten, averaged over trajectory. 
J(E =145 Mev) =16 +3 protons/cm’ sr sec 





Integral omnidirectional proton flux at point based on satellite data.° 


























proximately equal to the number in the back- The spectra obtained at the three points are 
ground area in the same energy range. There- shown in Fig. 3. 
fore, it is not possible, with the statistics avail- Three things are apparent from these results: 
able at this time, to make reliable measurements (1) At the higher latitudes, the slope of the spec- 
of the spectrum above about 100 Mev. trum below 40 Mev is very steep compared to 
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predictions from galactic cosmic-ray neutron 
albedo theory.* (2) The shape of the spectrum 
and the flux change with latitude. Points 3 and 5 
are at the same altitude and geomagnetic longi- 
tude. (3) At a comparable position in the belt 
(point 5) the flux and shape of the spectrum agree 
with the previous data. 

Singer* has predicted a change in the slope of 
the spectrum at higher energies with a change 
in latitude based on the assumption that the max- 
imum energy of the trapped protons decreases 
with altitude due to the breakdown of the adia- 
batic invariance of the magnetic moment. It is 
also possible that the change in the spectrum at 
lower energies is due either to solar protons 
trapped in the geomagnetic field or to protons 
produced by the decay of albedo neutrons pro- 
duced in nuclear interactions of solar protons 
over the polar cap. In either case there would 
be a preferential trapping of low-energy particles 
at the higher latitudes. 

It is not clear which of these mechanisms is 
responsible for the change in the spectrum. Fur- 
ther analysis will be made in these emulsions to 
determine the maximum energy at which protons 
can be trapped at a given point. Additional flights 
are scheduled in late 1961 to determine if thereis 
a time variation. Assuming that the low-energy 
protons are associated with solar activity, their 
flux should decrease as the number of solar pro- 
ton events per year decreases during the solar 
cycle. Changes in the spectrum after a large 
flare have been reported by Armstrong et al.? 

The unidirectional integral flux, J, for protons 
with E 231 Mev is given in Table I. In order to 
compare this flux with the measurements in sat- 
ellites, the assumption has been made that the 
flux is isotropic so that the integral omnidirec- 
tional flux is obtained by multiplying J by 47. 





The values obtained in this way are listed in 
Table I together with the best value of the flux 
at the three points based on the analysis of the 
satellite data.* The quoted errors are large due 
to the effects of the background on the integral 
flux. The agreement is quite good considering 
the errors. The unidirectional flux, assuming 
isotropy, is consistently higher than the omni- 
directional flux at all three points, indicating 
that the assumption is not valid. 

The nose cone was spin stabilized nearly paral- 
lel to the magnetic field lines, and the unidirec- 
tional flux in the emulsion was measured pri- 
marily at angles near 90° with respect to the 
field lines. At low altitudes, where there is a 
high density of mirror points, the flux should 
not be isotropic but should be higher at 90° an- 
gles as observed. 

The solution of the problems associated with 
the exposure and recovery of emuisions at 1880 
km required outstanding and unique contributions 
from a number of people. These efforts are 
gratefully acknowledged by the authors. In par- 
ticular, the shot would not have been successful 
without the work over and above his normal re- 
sponsibility by Project Engineer Charles E. 
Campbell, of the Goddard Space Flight Center. 





1§. C. Freden and R. S. White, Phys. Rev. Letters 
3, 9 (1959). 

2a. H. Armstrong, F. B. Harrison, H. H. Heckman, 
and L. Rosen, J. Geophys. Research 66, 351 (1961). 

35. C. Freden and R. S. White, J. Geophys. Research 
65, 1399 (1960); S. F. Singer, Phys. Rev. Letters 1, 
141, 181 (1958); W. N. Hess, Phys. Rev. Letters 3, 
11 (1959). 

4s. F. Singer, Space Research (Interscience Pub- 
lishers, Inc., New York, 1960), p. 797. 

5C. E. McIlwain and J. A. Van Allen (private com- 
munication). 
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FERROMAGNETISM IN SOLID SOLUTIONS OF SCANDIUM AND INDIUM 


B. T. Matthias, A. M. Clogston, H. J. Williams, E. Corenzwit, and R. C. Sherwood 
Bell Telephone Laboratories, Murray Hill, New Jersey 


(Received June 12, 1961) 


We have found that a solid solution of approxi- 
mately 24 atomic percent indium in scandium 
becomes ferromagnetic below 6°K. The ferro- 
magnetism occurs over an extremely narrow 
range of compositions between Sc, In, 25, and 
SCo.rselMo.242, & Change of 0.4 atomic percent of 
indium. 

The ferromagnetism of these solid solutions 
was encountered during a search for supercon- 
ductivity in the scandium-indium system. After 
the discovery that La,In becomes superconduct - 
ing at 10.4°K,' the corresponding phases of yttri- 
um and scandium became of interest. While In 
is not soluble in Y, no difficulty has been found 
in dissolving In in Sc simply by melting the ma- 
terials together in an arc furnace. Measure- 
ments of the susceptibility of these melts in- 
itially showed a maximum near 20 atomic per- 
cent In. After the ingots were annealed for a 
week at 800°C the maximum susceptibility in- 
creased greatly and moved to higher In concen- 
trations. No superconductivity has been found 
in this system. X-ray patterns do not show the 
cubic La,In phase but instead the essentially 
hexagonal close-packed symmetry of scandium 
metal, with a smaller lattice constant, consist- 
ent with a smaller atomic radium for metallic 
indium than for scandium, At the composition 
SCo7olNMp,39 the structure changes. 

In Fig. 1 we show a plot of the magnetic mo- 
ment in cgs units per gram observed in a field 
of 14000 gauss at 1.4°K as a function of the 
atomic fraction of In for the annealed samples. 
The extreme sharpness of the peak shows how 
very easy it would be to miss the occurrence 
of ferromagnetism in this system. Measure- 
ments have been made of the susceptibilities x 
for these solid solutions from 1.4°K to room 
temperature. x is strongly temperature de- 
pendent, becoming large at low temperatures, 
and can be fitted approximately to a Curie-Weiss 
law of the form x=Nup*pese*/3k(T-0) where N 
is the number of magnetic centers, up is the 
Bohr magneton, peff is the effective magnetic 
moment per atom, and @ is the Curie tempera- 
ture. The values of Pegg and 6 are shown as a 
function of composition in Figs. 2 and 3 and in 
Table I. Although por¢ Shows a small peak at the 
critical composition, the large peak in the mag- 
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FIG. 1. Magnetic moment per gram observed at 
1.4°K in a field of 1400 gauss with a sample of 
Scj-,In, as a function of atomic fraction of In. 
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FIG. 2. Effective magnetic moment in Bohr magne- 
tons per atom for Sc j-xIny as a function of atomic 
fraction of In, determined by fitting observed sus- 
ceptibilities to a Curie-Weiss law. 
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FIG, 3. Curie-Weiss constant @ in degrees Kelvin of 


of a sample of Sc,_,In,, as a function of atomic frac- 
tion of In, determined by fitting observed suscepti- 
bilities to a Curie-Weiss law. 


netic moment shown in Fig. 1 clearly arises from 
the fact that @ becomes positive at this point. The 
accuracy in determining @ is not better than 1 or 
2°K. 

A positive value of 6 implies that the material 
must become ferromagnetic for T< 6. For the 
four samples where @ is positive, the magnetic 
moment has been studied as a function of H and 
T. The Curie temperature Tc; has been deter- 


mined by plotting H/o, vs 0? isotherms, where Og 


is the magnetization ler gram. These results are 


given in Table II along with the saturation mo- 
ment ug per atom in Bohr magnetons. We also 
show the magnetic moment per atom determined 
from the Curie-Weiss law assuming g=2. Rema- 
nence has been observed at 1.4°K. 

The nature of the ferromagnetic state in 
SCo,76IMo,24 is not clear at this point. The Curie- 
Weiss law obeyed by the susceptibilities suggests 
that the scandium d electrons are partially lo- 
calized and exhibit a small magnetic moment. 


Table I. @ and poss for Sc, _,In, as a function of x. 








a 4 

x beff aK 
0.222 0.60 -1l 
0.235 0.60 -10 
0.238 0.66 +5 
0.240 0.70 +6 
0.241 0.69 +6 
0.242 0.65 +7 
0.244 0.62 -6 
0.250 0.56 -5 
0.267 0.50 -15 











If this were the case we would have a ready ex- 
planation of the fact that no superconductivity 
develops in this system in contrast to the high 
superconducting transition temperature of La, In. 
This model, however, provides no explanation 
why ferromagnetism should appear in such an 
extremely narrow range of compositions. 

The sensitivity to composition and homogeneity 
is perhaps suggestive of a band-filling mechanism, 
leading to itinerant electron ferromagnetism. A 
similar situation has been reported previously 
for ZrZn,,? which becomes ferromagnetic below 
35°K. 

In this case, Abrahams*® could not detect the 
1.4 Bohr magnetons per atom indicated by the 
Curie-Weiss law although his measurements 
were consistent with the smaller value of 0.17 1p 
determined by measurements of the saturation 
magnetization. Electrical resistivity measure- 
ments by Olsen‘ also failed to show any spin dis- 
order scattering above the Curie temperature. 
The ferromagnetism of ZrZn, is easily destroyed 
by impurities and other lattice disturbances. 

An attempt to account for the ferromagnetism 
of the Sc-In system by an itinerant electron mod- 
el, however, also faces difficulties. It is easy to 
show that the susceptibility of the band model 
above the Curie temperature will be given, in the 
Hartree-Fock approximation, by the expression 


4.7... 2 J, (1) 


where x,(7) is the temperature -dependent suscep- 
tibility that would prevail in the absence of ex- 
change, 2 is the atomic volume, andJd is approxi- 
mately the exchange self-energy of an electron on 
one lattice site. If the band in question is fairly 
wide, x,(7) and x will be essentially temperature 
independent, contrary to our observations. If the 
band is narrow enough to introduce substantial 
temperature variations in x,(7), it will still be 
difficult to match the shape of the observed sus- 
ceptibility curves using Eq. (1) since experimental- 


Table Il. Tc and pg for Sci; ~~)In, as a function of x. 





x 1 ,°K 4y,(Bohr magnetons) (1+ P o¢e?)?-1 





Cc 
0.238 5 0.035 0.20 
0.240 6 0. 036 0.22 
0.241 6 0.039 Q.22 
0.242 6 0.038 0.18 
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ly x is rising sharply at the lowest temperatures, 
The narrow range of compositions for which 
ferromagnetism exists in the Sc-In system sug- 
gests that this kind of ferromagnetism may be a 
more general phenomenon than has been thought 
until now, and may be found to occur in other 
systems. 
We would like to thank Mrs. V. B. Compton 





for her x-ray investigations of these alloys. 
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SHIFT OF NUCLEAR QUADRUPOLE RESONANCE FREQUENCY BY ELECTRIC FIELD 


Toshimoto Kushida* and Kuniomi Saikit 


Hiroshima University, Hiroshima, Japan 
(Received June 14, 1961) 


It is well known that a static magnetic field 
splits the nuclear quadrupole resonance (NQR) 
lines in a single crystal.’ One might immediately 
think of the similar effect caused by an electric 
field applied to the specimen. As a matter of 
fact, Pound? applied a field of about 50000 v/cm 
across a thin plated-shaped single crystal of po- 
tassium iodide, but no shift in the resonance fre- 
quency nor change in the line shape was observed 
in the I’*” nuclear resonance line. Gutowsky and 
Williams® tried to observe the effect of an elec- 
tric field on an NQR line of the Cl nucleus in 
NaClO, single crystal without success. 

Recently Bloembergen* pointed out the possibil- 
ity of observing the change in the quadrupole 
coupling constant caused by an external electric 
field, if the crystal lattice has no local inversion 
symmetry about the pertinent nucleus. 

It is the purpose of this Letter to report the suc- 
cessful observation of this effect for the Br®™ NQR 
line in a single crystal of NaBrO,.5 The resonance 
frequency and the line shape are observed as a 
function of applied electric field using a regen- 
erative spectrometer® in conjunction with a pre- 
cise 0°C ice bath in order to keep the tempera- 
ture of the sample constant during the run. Alu- 
minum film electrodes are evaporated on the sur- 
faces of a thin plate sample which has been cut 
from a larger single crystal. 

The peak frequency of the line shifts linearly 
with the applied field E when E£ is parallel to 
[111], whereas no shift in the center frequency 
is observed in the case of £||[100]. Reversal of 
the polarity of E reverses the sign of the shift. 

A typical example of the observed absorption 
curve is shown in Fig. 1. Here E is along the 
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FIG. 1. The effect of the electric field on the Br*! 
resonance in NaBrO; single crystal. The arrow 4 
indicates the center frequency of the original line 
without the field. 


[111] direction of the crystal, and H,, an rf 
magnetic field in a sample coil of the spec- 
trometer, is applied along the [110] direction. 
The vertical arrow ; indicates the center fre- 
quency without the electric field. The line shape, 
which was symmetric without an E£ field, is slight- 
ly asymmetric, and the center frequency has 
shifted by about 500 cps. 

A cubic unit cell’ of NaBrO, has four physically 
inequivalent, but chemically equivalent, BrO,~ 
radicals. The maximum principal axes of elec- 
tric field gradients about Br nuclei are parallel 
to the body diagonals of the unit cell. The field 
gradients along the [111], [111], [111], and [171] 
directions will be denoted, respectively, by q,, 
Ja» Ys, and q,. 

The distortions of the field gradients produced 
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by the electric field along [111] are the same for 
Go, Ys, and q, because of the threefold symmetry 
of the crystal along this direction. Because q, 
distorts differently, a doublet separation of the 
line would be expected. Since the relative inten- 
sity of the component lines is proportional to 
sin* 2(H,,q), the intensity ratio of the compo- 
nents of the doublet will be /,:/,,,=3:5 in the case 
of H,||[110]. Using this intensity ratio the ob- 
served asymmetric curve can be decomposed 
into two components as shown in Fig. 1. When 
H, is along the [111] direction, one observes 
only the (¢,9,9,) line, which does not show any 
asymmetry nor increase in the linewidth com- 
pared with the original line without the E field. 
The observed shift of each component is 


q, line -2500+ 500 cps/(10 000 v/cm), 
(429394) line 750+ 150 cps/(10000 v/cm). 


The minus sign for the q, line indicates only the 
fact that the q, and the (g,q,9,) lines shift to op- 
posite directions; no attempt has made to de- 
termine the absolute sign of these magnitudes. 

When E is along [100] and H, is along [010], 
the shift for the (q,q,) line is the same as, but 
opposite in sign to, that for the (q,q,) line. Since 
both lines have the same intensity, no shift in 
the center frequency is observed up to about 
5000 v/cm, where the crystal always breaks 
down. 

The distortion of the field gradient tensor com- 
ponent, OP 57» may be expressed as 


The third rank tensors, 89;;/8E p, are zero if 
local inversion symmetry about the nuclei exists. 
This is the case for KI. Symmetry requirements 
of the crystal structure® and the traceless prop- 
erty’ of ij generally decrease the number of the 
independent components of this tensor. 

The observed shift of the resonance frequency 
indicates only the change in g, which is 9;; along 
the maximum principal axis, because a small 
asymmetry parameter 7 produced by the E field 
changes the resonance frequency only in the sec- 
ond order. From the shift of the q, line we ob- 
tain 

a /Vo =(4 “*9q/8E )E, (2) 


and the shift of the (q,q,q,) line gives 


ava x . 
6 = “lp E Breda " 1 z 
Y anqeag 07 49 24/8E gE ~ (4g °84/8E y)3 


(3) 
10 


where vy, and q, are, respectively, the original 
resonance frequency and q without the field. 

The values of the “field-gradient electrostatic” 
tensor are derived as 


qo 28q/8E “7 (1.7+ 0.3) x1075/(10 000 v/cm), (4) 


and 


o '0q/8E | = (0.04 0.2) x1075/(10 000 v/em). (5) 


The vanishing of 8q/aE, is expected from the sym- 
metry requirements of the third rank tensor for 
this crystal.*”® 

Since NaBrO, is piezoelectric,® it would be in- 
formative to compare the value of g,~*8q/2E , 
with the q change produced by a piezoelectric 
strain. The electric field along the [111] direc- 
tion produces an axial strain along this direction.® 
The resonance frequency change accompanying 
this axial strain can be estimated from the piezo- 
electric constant and the pressure (hydrostatic) 
dependence of the resonance frequency.’® Since 
the implicit effect via the piezoelectric strain 
thus estimated is less than 10% of the observed 
value of g,-'8q/8E |, it may be concluded that 
the explicit electric effect on g, namely 
(q.*8q/8E |) const strain’ 48 Mainly responsible 
for the observed effect. 

The authors are indebted to Professor N. 
Bloembergen for valuable discussions and to 
Dr. A. H. Silver for the critical reading of the 
manuscript. 
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LINEAR EFFECT OF APPLIED ELECTRIC FIELD ON NUCLEAR QUADRUPOLE RESONANCE* 


J. Armstrong, N. Bloembergen, and D. Gill 
Gordon McKay Laboratory, Harvard University, Cambridge, Massachusetts 
(Received June 14, 1961) 


Experimental evidence for the existence of a 
linear Stark effect which was predicted theo- 
retically by one of us,’ has been found in the 
zero-field quadrupole resonance of Cl®* in KC10, 
and NaClO,. Kushida? has independently, and 
somewhat earlier, discovered the same effect 
for the Br®™ resonance in NaBrO,. 

Because single crystals were not immediately 
available, our experiments have been carried out 
with powders of fine crystallites, which filled 
the 1.68-mm wide spacings between a sandwich 
of condenser plates. The tank coil of a modified 
Pound spectrometer, consisting of a few turns 
of copper ribbon, was wound around the electro- 
static condenser assembly and placed in liquid 
nitrogen. A dc voltage of up to 3500 volts could 
be applied without causing electrical breakdown. 
This corresponds to an average field strength in 
the crystallites of about 17.1 kv/cm. This value 
was determined from approximate formulas for 
the inhomogeneous medium‘ with the known pow- 
der filling factor of 0.63. 

The frequency of the spectrometer was swept 





FIG. 1. The zero-field 
quadrupole resonance of Cl® 
in polycrystalline KC1O; at 
77°K with and without applied 
electric field. 





linearly through the resonance, which occurs 

at 28954 Mc/sec for KC1O, at 77°K. A frequency 
modulation at 340 cps with a modulation index 

less than 0.25 was used.® The signal output of 

the phase-sensitive detector corresponds to the 
absorption derivative of the zero-field quad- 
rupole resonance. Two recorded traces for KC1O,, 
one without electric field and one with 3000 volts 
between the condenser plates, are shown in Fig. 1. 
The marked variation in line shape is striking. 
From such reproducible recordings the increase 
in second moments as a function of the applied 
electric field was calculated. The results in 
Table I show that the second moment increases 
linearly with the electrical energy stored in the 
crystallites. 

The theoretical interpretation of these results 
is based on a linear variation of the electric 
field gradient at the position of the nuclei with 
the externally applied electric field, if these 
positions lack inversion symmetry. This rela- 
tionship can be described by a third order ten- 
sor with components Rjjz,. The variation in field 
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Table I. The observed second moment (A:*) and the 
width 6v between points of extreme slope of the cis 
resonance as a function of the electric field strength E 
in polycrystalline KC1O, and NaClO;. 

















Einkv/em (Av) inecps? _ dv in cps 
KCI1O; 0 65 000 480 
7.32 91000 530 
14.6 137000 990 
NaClO, 0 105 000 610 
7.32 117000 690 
12.2 131000 785 
17.1 eee 905 
18.3 151000 eee 
gradient components is then given by 
= Oe Baa Os 
e(aq)= Le Ry Ept LL Si etm 
R=X,y,2 L,m=X,y,2 
(1) 


The effect of variation with strain described by 

a fourth order tensor S was introduced earlier.*:’ 
In piezoelectric crystals the application of an 
electric field produces a strain 


Com pe etme R @) 
The interest here is clearly in the new physical 
effect of an electric field at constant strain. The 
symmetry properties of third and fourth order 
tensors are well known.® In fact, the relationship 
given in Eq. (1) is quite similar to the relations 
between the electro-optic and elasto-optic effects. 
In our case the symmetry at the location of the 
individual nucleus counts rather than the sym- 
metry of the crystalline unit cell. Since the elec- 
trostatic potential obeys Laplace’s law, some 
trace relations exist, 


UR, p= US 


ii 1 iilm 


0. (3) 


In the most general (triclinic) case R has fifteen 
independent constants. 

The crystal KC1O, has monoclinic symmetry 
with inversion center (Coy or 2/m). The crystal 
is not piezoelectric and has two equivalent mole- 
cules in the unit cell. The symmetry at the posi- 
tion of the Cl nuclei is that of a single mirror 
plane (C; or m). The electric field effect is in- 
deed observed at constant strain. In each crys- 
tallite the zero-field splittings of the two Cl nu- 
clei in the unit cell are displaced by equal and 
opposite amounts. The Kramers degeneracy of 


12 


the nuclear spin levels is not lifted. The effect 
has nothing to do with the existence of a nuclear 
dipole moment.’°*!_ From a single observation 
of the variation of zero-field splitting, the eight 
independent components of the R tensor in the 
Cs symmetry can clearly not be determined. 
With a single crystal one could measure the 
variation in the zero-field splitting with the 
applied electric field in two directions with- 

in the plane of symmetry. This would still give 
only two components. The other components 
can in principle be determined by two methods 
which will be described below. 

The crystal symmetry of NaClO, is cubic (23 
or T). The symmetry at the position of the four 
Cl nuclei in the unit cell is trigonal (3 or C;,). 
The crystal is piezoelectric. It can, however, 
be estimated from the known piezoelectric con- 
stant and the known effect of hydrostatic com- 
pression on the quadrupole splitting’* that the 
observed effect is mostly due to the first term 
on the right-hand side of Eq. (1). The first order 
change in the zero-field splitting v, for this case 
of axial symmetry is thus given by 


Av/v9= (49, ,/19) =Ro22” 2/9 (4) 
where the z axis is taken at each Cl nucleus 
pointing to the nearest Na nucleus. The increase 
in second moment in the powder is 


~ Eo 9/49" (5) 


oil 
(4 ar R22 


The same formula may be used in a good approxi- 
mation for KC10,, since the deviations from axi- 
al symmetry in that case are quite small. This 
is borne out by the small value of the asymmetry 
parameter,’® 7=0.007. There can be little doubt 
that a single crystal of KC1O, would show a re- 
solved doublet for the two Cl nuclei, which be- 
come nonequivalent on application of an electric 
field. This would throw further light on the 
structure of this crystal.'* 

The value of |Rz2z21 and |A4q,,| derived from 
experiment with Eqs. (4) and (5) agrees well with 
a simple theoretical estimate based on a point 
charge model. The electric field displaces the 
nearest alkali ion with respect to the Cl nucleus. 
The displacement Ar can be estimated either 
from the force constants of the ionic crystals 
or from the contribution of ionic displacements 
to the static dielectric constant € 4, 


= — 
Ar (€4. n®*)E /4N@ | o> 
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where N is the number of positive ions per unit 
volume with effective charge eeg¢ and n is the 

optical index of refraction. The change in quad- 
rupole coupling is in order of magnitude given by 


e4q = 6(1 “Va ogg” 47> 


or 
R,,, 78 -v, eq, -m4mNYA), 6) 

where 1 -y,, is the Sternheimer factor, andr 

the anion-cation distance. 

This oversimplified model should, of course, 
be improved in several ways, but it gives the 
observed value |Rz22 | =2.1 10° cm™ with 
1 -y.. =35. Its purpose is merely to illustrate 
the physical origin of the observed R tensor, 
which is closely related to the third derivatives 
of the electrostatic potential. 

It would appear that only one constant, R,,,, 
can be measured in the case of axial symmetry, 
both in the powder and in a single crystal. This 
would leave four other independent constants of 
the R tensor for C, symmetry in NaClO, unde- 
termined. They relate to off-diagonal elements 
in the electric quadrupole interaction. If, how- 
ever, an external magnetic field is applied to a 
single crystal at an angle with the axis of sym- 
metry of sufficient strength to scramble the nu- 
clear spin states so that the m,’s are poor quan- 
tum numbers, the various resonant frequencies 
between the four nuclear spin states will vary 
linearly with applied electric field components 
in the x and y direction. In principle, all com- 
ponents of R may be determined from such meas- 
urements. 

It is perhaps of more interest to utilize the ex- 
istence of off-diagonal elements proportional to 
the applied electric field for the induction of nu- 
clear spin transitions by a radio-frequency elec - 
tric field. The transition probabilities are pro- 








portional to the square of the electric field 
strength and can readily be calculated from Eq. 
(1) and the known matrix elements of the nuclear 
quadrupole operator. The components 4q,, and 
Aqyz will be connected with Am ,=+1 transitions 
and Aq,y -Adyy with Am ;=+2 transitions. The 
resonant absorption may be described by an ef- 
fective nuclear dielectric susceptibility, «”/47. 
The coupling strengths of the nuclear spin to ap- 
plied electric and magnetic fields may be com- 
pared with each other. The last column in Table 
II has only order-of-magnitude significance since 
differences in magnitude between the nonvanish- 


ing components of R and between the nuclear 
spin operator matrix elements have been ignored. 


In principle, accurate measurements of the ef- 
fective nuclear dielectric susceptibility tensors 
would determine the other components of R. 
Such experiments bear, of course, a close re- 
semblance to acoustically induced resonances.’ 
The electrical resonances at constant strain 
have the advantage that acoustical propagation 
effects are absent. Electric radio-frequency 
fields of the required magnitude and uniformity 
can readily be obtained. 

For the heavier elements the coupling to the 
electric field becomes large enough that it may 
be detectable directly as a decrease in Q ofa 
capacitance. An electric field of 3 kv/cm ap- 
plied at the Cl°* resonant frequency in KC1O, at 
room temperature would saturate that spin reso- 
nance. Double-resonance experiments are feasi- 
ble, in which, e.g., an electric radio-frequency 
field saturates a |Am|=2 transition and a mag- 
netic radio-frequency field monitors a |Am|=1 
transition. Electric 90° and 180° pulses are 
possible. Numerous other electric variations 
of magnetic resonance experiments suggest 
themse!ve 3. 

The phenomenological description given here 


Table II. Comparison of electric and magnetic nuclear spin coupling. 








Quadrupole Stark shift for 104 v/em 
splitting along axis of symmetry 
Nucleus Crystal vo in Mc/sec Avineps Ad in gauss €’’ /4ny’’ 
0 
Na’ NaClO, 0.89 50° 0.044" 2x10-6# 
cis NaClO; 30.6 250° 0.6 3.3x10~44 
Br*! NaBrO; 150 2500° 0.23 5 x1073* 








*-rheoretical estimate. 
From data of Table I and Eqs. (4) and (5). 
Experimental value from reference 2. 
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is quite general and can readily be applied to 
other situations. Consider, e.g., integral spins 
in axial symmetry. The levels m, and -m, are 
degenerate. An applied electric field off axis 
will lift the degeneracy of the levels m,=+1 in 
first order. This is a true linear Stark splitting 
of a degenerate level. 

An interesting situation occurs if the nucleus 
occupies a cubic site without inversion symmetry. 
In the III-V compounds with ZnS structure the 
nucleus occupies a site with Tg= 4 3m symmetry. 
There is no quadrupole interaction and an exter- 
nal magnetic field will give 2] +1 equally spaced 
levels with one resonant frequency. If an electric 
field is applied— some III-V compounds, such as 
GaAs, can be obtained of sufficiently high resis- 
tivity that fields of 10 kv/cm or more can be 
maintained—the nuclei are displaced from their 
positions of cubic symmetry. A first order quad- 
rupole splitting or broadening proportional to E 
should result, which can readily be calculated 
for arbitrary orientations of E and H,. The R 
tensor is characterized by one constant in this 
case. 

Similar effects in electron spin resonance can 
be expected, if a paramagnetic ion occupies a 
site without inversion symmetry. Odd terms in 
the crystalline field potential will give rise to 
changes in the constants of the spin Hamiltonian 
on displacement of the ion by an applied electric 
field. Linear variations in optical absorption 
frequencies of transition metal ions may also 





occur.’* The ratio of the displacement of the 

resonant line for a given applied electric field 
to the linewidth is generally less favorable in 

these cases. 
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DIRECT EXPERIMENTAL MEASUREMENT OF THE MAGNETIC FIELD DEPENDENCE 
OF THE SUPERCONDUCTING ENERGY GAP OF ALUMINUM 
D. H. Douglass, Jr. 


Lincoln Laboratory,* Massachusetts Institute of Technology, Lexington, Massachusetts 
(Received June 7, 1961) 


This Letter reports experimental measurements 
of the field dependence of the superconducting en- 
ergy gap of aluminum. The energy gap was meas- 
ured directly by using the electron tunnelling tech- 
nique.’*? The magnetic field was longitudinal and 
the reference superconductor was lead. In Fig. 1 
the normalized energy gap vs reduced magnetic 
field is shown for two films of thickness 3000 A 
and 4000 A. It should be noticed that for the 3000A 
film the energy gap goes smoothly to zero as the 
critical field is approached. It is described quite 
well by the equation [€(H)/e(0)]*=1-(H/H,)*. This 
behavior implies a second order phase transition. 


14 


Measurements on two other films of thickness 
500 A and 2000 A gave similar results. Also, 
these results were independent of temperature 
over the range available (0.75<7/7,.<1.00). For 
the 4000A film an entirely different behavior is 
observed. The gap initially drops less rapidly 
with field than in the previous case, and as the 
critical field is approached the energy gap drops 
very abruptly from a finite value to zero; this 
abrupt change usually occurs within less than 1% 
of the critical field. The value of the gap just 
prior to this abrupt drop will be called the crit- 
ical gap. The critical gap of the 4000A film was 
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FIG, 2. Energy gap vs magnetic field according to 


FIG. 1. Energy gap.of aluminum vs magnetic field 
for films of thickness 3000 A and 4000 A. Two other 
films of thickness 500 A and 2000 A behaved similarly 
to the 3000A film. The dashed curves are the best 
straight lines through the data points. 


found to be a function of temperature. 

Douglass® has pointed out that the Ginzburg- 
Landau* theory with Gor’kov’s® addition could 
be used to explain the field dependence of the 
energy gap. Figure 2 shows the predictions of 
the theory which were obtained by solving Eqs. 
(1) and (2) of reference 3. According to the the- 
ory, the nature of the curve is a function only 
of the ratio of the thickness, d, to the penetra- 
tion depth, A. For small values of d/d the curves 
obey the equation [€(H)/e(0)]?=1-(H/H,)*. As 
d/x becomes larger the curves bow upward with 
the energy gap continuing to go smoothly to zero 
at the critical field (i.e., the critical gap remains 
equal to zero). For d/x>V5 the critical gap as- 
sumes a finite value, indicating that the phase 
transition has become first order. Thus, for d/A 
<V5 the theory predicts a critical gap of zero and 
a second order phase transition, and for d/A > v5 
it predicts a finite critical gap and a first order 
phase transition. The experimental data of Fig. 1 
show that these two types of predicted behavior 
are in fact observed, although the predicted bow- 
ing of the curves does not seem to occur. Values 
of d/a for the various films were calculated® from 
the measured critical fields using the Ginzburg- 
Landau theory. For a given thickness one can 


the Ginzburg- Landau-Gor’kov theory for various ratios 
of thickness to penetration depth. 


vary d/d over a range of values by changing the 
temperature since A is temperature dependent. 
The largest attainable value of d/a is limited by 
the value of A at low temperatures, and the small- 
est value of d/A is determined by how close one 
can get to the critical temperature and still have 
signals big enough to measure. Values of d/A be- 
tween 2.42 and 3.62 for the 4000A film were ob- 
tained by varying the reduced temperature between 
0.749 and 0.965. 

Experimental measurements of the critical gap 
are shown along with the theoretical curve’® in 
Fig. 3. Measurements were made on four dif- 
ferent films covering a range of d/A from 0.10 
to 3.52. It is seen that the data agree very well 
with the theory when d/a is less than 2.8. The 
transition from first to second order is not in- 
consistent with d/A =V5; the data show that it 
occurs for a d/d between 1.9 and 2.4. For d/A 
greater than 2.8 the value of the critical gap de- 
viates from the theoretical curve. 

There are two limitations on the Ginzburg- 
Landau-Gor’kov theory. The first is that it is 
true only in the local limit. This means that the 
penetration depth A must be larger than the co- 
herence length ~, which certainly is not true in 
the bulk state of aluminum where A(T =0)=500 A 
and ,=15000 A. However, one would expect 
that for a thin evaporated film the coherence 
length would be limited by the size of the crys- 
tallites. Thus the condition of locality would be 
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4 >d, which is much easier to satisfy, and non- 
local effects will start to become important when 
d becomes comparable to A. Therefore, it is not 
surprising that deviations from the theory occur 
at d/A=3. For a metal for which the bulk pene- 
tration depth is greater than the bulk coherence 
length, the theory should be good for all d/A. 
The second limitation concerns the fact that 
Gor’kov derived the Ginzburg-Landau equations 
from the microscopic theory under the condition 
that the energy gap was small compared to its 
value at zero temperature and zero field. The 
theory is therefore valid near the transition tem- 
perature where the gap is small. For thin films 
it can be made valid even at T =0 by turning ona 
magnetic field. 

In summary, we have measured the magnetic 
field dependence of the superconducting energy 
gap of aluminum films as a function of the ratio 
of the film thickness to the penetration depth. 

The critical gap was found to be zero for values 
of this ratio less than 1.9 and equal to a finite 
value for values greater than 2.4. The agreement 
between these experiments and the Ginzburg- 
Landau-Gor’kov theory is very good. 

The author is grateful to R. Tweed for the fab- 
rication of the films and to Dr. R. Meservey, who 
assisted in some of the data taking. He is also 
grateful for helpful discussions with Dr. E. Max- 
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well and Dr. C. E. Chase. 


Note added in proof. Morris and Tinkham,’ 
from thermal conductivity measurements, have 
also determined the field dependence of the en- 
ergy gap in thin films. They come to conclusions 
similar to those presented here. 








*Operated with support from the U. S. Army, Navy, 
and Air Force, 

‘I, Giaever, Phys. Rev. Letters 5, 464 (1960). 

J, Nicol, S. Shapiro, and P. H. Smith, Phys. Rev. 
Letters 5, 461 (1960). 

3p. H. Douglass, Jr., Phys. Rev. Letters 6, 346 
(1961). 


‘yv, L. Ginzburg and L. D. Landau, J. Exptl. Theoret. 


Phys. (U.S.S.R.) 20, 1064 (1950). 

°L. P. Gor’kov, J. Exptl. Theoret. Phys. (U.S.S.R.) 
36, 1918 (1959) (translation: Soviet Phys.—JETP 9, 
1364 (1959)]. 

®There are as yet no data relating d/d tod for alu- 
minum. It might be argued that using the Ginzburg- 
Landau theory to calculate d/A from the critical field 
amounts only to a self-consistency check; however, 
estimates of d/A using the approximate thickness and 
the bulk penetration depth, modified by size effects ac- 
cording to M. Tinkham [Phys. Rev. 110, 26 (1958)], 
give values which do not differ by more than 20-30% 
from those quoted here. 

'D. E. Morris and M. Tinkham, Phys. Rev. Letters 
6, 600 (1961). 
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PARAMAGNETIC RESONANCE DETECTION OF LUMINESCENT CENTERS AND TRAPS 
IN SELF-ACTIVATED ZnS PHOSPHORS 


Paul H. Kasai and Yoshiro Otomo 


Hitachi Central Research Laboratory, Kokubunji, Tokyo, Japan 
(Received May 22, 1961) 


The atomic nature of the luminescent center 
in self-activated ZnS phosphors has been the 
subject of many papers.** Kréger and Vink 
have proposed that the self-activated centers 
are Zn** ion vacancies that have lost one elec- 
tron from their nearest surrounding S~~ ions. 
Magnetic susceptibility measurements by Bowers 
and Melamed,’ however, revealed that the self- 
activated center is diamagnetic, whereas the 
model proposed by Kréger and Vink should be 
paramagnetic. Prener and Williams® then sug- 
gested a model in which the centers in question 
are Zn** ion vacancies with one of their charge- 
compensating ions (Cl~, Br~, I~, Al***, or Ga***) 
associated at the nearest possible sites in the 
ZnS lattice. More recently Prener and Weil* 
supported the latter model by showing experi- 
mentally that the emission peak of the self-acti- 
vated luminescence depends on whether a group 
III-B or group VII-B element is used as a coac- 
tivator. The existence of various traps for ex- 
cited electrons or holes in luminescent ZnS phos- 
phors including the self-activated type are also 
well known,°® though here again the problems con- 
cerning the exact nature of some of these traps 
are still unsettled. If the self-activated ZnS 
phosphor is truly diamagnetic in its unexcited 
state, then the presence of such traps implies 
that the phosphor should become paramagnetic 
when enough electrons are transferred from the 
luminescent centers into these traps as is often 
done in thermoluminescence experiments. 

The following is the preliminary report of our 

electron paramagnetic resonance (EPR) detec- 
tion of such excited states frozen at the liquid 

nitrogen temperature. 

The polycrystalline samples of self-activated 
ZnS phosphors were prepared by firing the lu- 
minescent grade ZnS powder using various al- 
kali halides (NaCl, NaBr, KCl, and KBr) as 
flux at 800°C for 3 hours. A comparatively low 
firing temperature was employed in order to en- 
sure the formation of cubic crystals. The x-ray 
analysis did prove that they were all cubic, and 
all of them exhibit a single emission band peak- 
ing around 4640 A characteristic of the self- 
activated emission. 

In parallel with the result of Bowers and Mel- 


amed, none of these samples give an EPR signal 
in the dark, either at room temperature or at 
77°K, except six equally spaced weak lines due 
to a trace of Mn** impurity. At room tempera- 
ture, simultaneous illumination of the sample 
with 3650A ultraviolet radiation produced no 
change in the observed spectra. At liquid ni- 
trogen temperature, however, the simultaneous 
ultraviolet illumination gave rise to two addi- 
tional signals A and B in every sample. Figure 
1 shows such spectra obtained with ZnS(NaCl) 
and ZnS(NaBr), respectively. Changing the ca- 
tion of the flux from Na* to K* has no effect upon 
the spectra. 

The shape of the signal B is more or less sym- 
metric, and is little affected when the halide 
ion in the flux is changed. The g value of the 
signal B is 1.883+0.001. The fact that its g 
value is less than 2 implies that the signal B 
is associated with un-ionized donors, namely 
filled traps.® 

The shape of the signal A is just what one would 
expect for a polycrystalline sample having an an- 
isotropic g value.” Also, as can be seen from the 
figure, when Cl” in the flux is replaced by Br’, 
the apparent linewidth becomes much broader. 
In either case the g values observed on the sig- 
nalA are 


g = 2.0524 0.001, 


&y = 2.027 + 0.005. 







ZnS (NaC!) 


ZnS (NaBr) 





FIG, 1. EPR spectra of self-activated ZnS phosphors - 
under ultraviolet illumination at 77°K. Six equally 
spaced weak signals are due to Mn** ion impurity, 
each spacing being approximately 68 gauss. 
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It should be noted that the mean g value is greater 
than the free spin value, implying that the signal 
A is caused by a hole, an empty acceptor. 

Both signals A and B are found to decay when 
the illumination is stopped. The rate of the de- 
cay is quite slow at 77°K except for the first few 
minutes immediately after stopping the excitation. 
Approximately half of the original intensity can be 
still detected after 5 hours, so long as the sam- 
ple is kept in the dark at 77°K. The rapid initial 
decay may be what is responsible for the glow 
peak observed near liquid nitrogen temperature 
in the thermoluminescence experiments.****® 
The subsequent slow but steady decay is probably 
due to the tunneling process discussed by Hoogen- 
straaten.°® 

Among many defects which can act as electron 
traps in the ZnS phosphors, the Co** ion is known 
to provide one of the deepest traps.° Consequently 
if our previously stated analysis concerning the 
nature of the signals A and B is correct, the 
introduction of Co** ions into these phosphors 
should result in a strong enhancement of the sig- 
nal A and weakening of the signal B. This was 
indeed found to be the case. Introduction of 10~° 
(mole ratio) CoCL, increased the intensity of the 
signal A by a factor of five while making the sig- 
nal B appear only as a trace. Figure 2 shows 
the spectrum obtained with ZnS(KC1) containing 
10-° Co**. A strong broad resonance at the left is 
due to Co** ions. The g value of this Co** reso- 
nance is 2.248+0.005.'° In Fig. 1 one also notices 
a weak negative dip near the g=2.00 position, the 
center of the six Mn** hyperfine components. 

This signal, however, cannot be incorporated 
into the signal A, since it does not grow upon 
the introduction of Co** ions. The nature of 
this signal is still unknown. 


A 


— H 


FIG. 2. EPR spectrum obtained with ZnS:10-5Co(KCl) 
under ultraviolet illumination at 77K. A strong broad 
resonance at the left is due to the Cott ion. Note the 
enhancement on signal A. The amplifier gain used for 
this record was about 1/3 of that used for the records 
in Fig. 1. 


18 






If the defect responsible for the signal A is a 
simple Zn** ion vacancy surrounded by four S~~ 
ions at the tetrahedral corners, the observed 
anisotropy in the g value and the dependency of 
the line shape upon the halide ion in the flux are 
both difficult to explain. We, therefore, believe 
that the signal A is caused by the defect proposed 
by Prener and Williams.* The Zn** ion vacancy 
surrounded by three S~~ ions and one halide ion 
at the tetrahedral corners could well account for 
the above-mentioned features of this signal. Also 
since there is no significant change in the g value 
when Cl in the flux is replaced by Br~, we must 
attribute the observed difference in the line shape 
to the difference in the nuclear magnetic mo- 
ments of the respective halogen nuclei. A careful 
study of the spectra revealed that the ratio of 
the apparent linewidth of the signal A (the width 
shown by the arrows in Fig. 1) of ZnS(NaCl) to 
that of ZnS(NaBr) is 1:2.8, in close agreement 
with the ratio of the average magnetic moments 
of the respective halogen nuclei. The nature of 
the defect responsible for the signal B is less 
certain. The chloride flux always produces more 
of this signal than the bromide, though there is 
no significant change in the line shape, nor in 
the g value. The sulfur vacancies produced by 
some chemical reaction involving the alkali hal- 
ide flux are probably responsible for this signal. 

The EPR spectrometer used for this work is a 
Varian 100-kc spectrometer, model V-4500, and 
the frequency of the microwave is about 9.2 kMc/ 
sec. 

A more detailed and fuller account of this work 
will be published in the future together with the 
work on a single crystal which is now in progress. 





‘Fr. A. Kroger and H, J. Vink, J. Chem. Phys. 22, 
250 (1954). 

*R. Bowers and N. T. Melamed, Phys. Rev. 99, 1781 
(1955). 

33. S. Prener and F. E. Williams, J. Chem. Phys. 
25, 361 (1956). 

43. S. Prener and D. J. Weil, J. Electrochem. Soc. 
106, 409 (1959). 

ow. Hoogenstraaten, Philips Research Repts. 13, 515 
(1958). . 

®See also J. Lambe and C. Kikuchi, J. Phys. Chem. 
Solids 8, 492 (1959). 

"J. W. Searl, R. C. Smith, and S. J. Wyard, Proc. 
Phys. Soc. (London) 74, 491 (1959). 

®N. Riehl and H. Ortmann, Z. Naturforsch. 10a, 896 
(1955). 


*°H. Arbell and A. Halperin, Phys. Rev. 117, 45 (1960). 


"See also F. S. Ham, G. W. Ludwig, G. D. Watkins, 
and H. H. Woodbury, Phys. Rev. Letters 5, 468 (1960). 
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MAGNETIC SHIELDING CONSTANTS IN HETEROPOLAR DIATOMIC MOLECULES 


C. William Kern and William N. Lipscomb 
Department of Chemistry, Harvard University, Cambridge, Massachusetts 
(Received June 5, 1961) 


The general theory of magnetic shielding of 
nuclei in molecules, first developed by Ramsey’ 
using perturbation theory, showed that the 
screening effect at the nucleus was due to mo- 
tions of molecular electrons. More recent vari- 
ational and atomic theories? of chemical shield- 
ing in simple molecules give results in agree - 
ment with Ramsey’s conclusions. As generally 
applied, the perturbation approach requires a 
knowledge of all the excited states of the mole- 





2n 


cule, including the continuum.* However, we 
show here that when a suitable choice is made 
for the vector potential, the excited state portion 
of Ramsey’s formula vanishes.* The magnetic 
shielding is then reduced to matrix elements 
over the ground-state wave function. 

When a diatomic molecule, AB, in a‘ ground- 
state configuration and containing 2m electrons 
is placed in a homogeneous magnetic field, H, 
the Hamiltonian for the electronic system can 
be written as® 


Yr 


1 yee 
KH =4 » ie 4 (tae, -080 432) +V, 


ual 


Au 





where R is the vector pointing from nucleus A to B, b is a dimensionless constant to be determined, 
and magnetic effects of nucleus B, molecular rotations, and electron spin have been neglected. 

Following Ramsey,’ we show that the expression for the magnetic shielding constant of nucleus A, 
averaged over all orientations of the molecule with respect to the applied field, becomes 


1 2n bR 2n 1 2 
ae oe 3 = - -1 
p=1 p=l1 n'=1 
| 2n 2n 2n 2n 
x1{0| >> M /y, %\n)Xn’| 2) M }0)+{0| >> M |\n’)<Xn'| >) M /r, *|0) 
uot xAp Ap unt XA UL pal XA nol xAp Ap 
2n 2n 2n 2n 

+ bR(O M 3in’)n’ 0) + bR(O n')<n' M ry, *10), 

{ fe teal ts In’)<n ponte ( Pre ie — xAu! je } 


where MyAy and py, are the x and y components of the angular and linear momentum of electron u. 
We have assumed right- and left-handed coordinate systems on nuclei A and B, respectively. Using 
the well known closure property,° the choice of 6 which renders the second-order perturbation sum 
zero is 








s 
fia (AE ,) ay 012s Meg Mead Aan 10) ms O1UM 4 Mear/A *10) 
i 4 ae J 
(AE). ROI2M A Pod” ie 10) ROILM, 4 Pool™ pe 10) 


where (AE})ay and (AE9) ay are average excitation energies corresponding to the pairs of operators 
(MyApu/YAp’»MxAp) and MyxAu/7A y*>Py yu)» respectively. We assumed (AE) ay=(AEQ)ay since, for 
heteropolar diatomic molecules, ali operators in the excited state sum belong to the same irreducible 


19 


VoLUME 7, NUMBER | 





PHYSICAL REVIEW LETTERS 





Jury 1, 1961 





Table I. Magnetic shielding in two heteronuclear diatomic molecules. 






a 














Calculated Experimental Calculated Experimental 

diamagnetic diamagnetic paramagnetic paramagnetic 

Molecule shielding of H shielding of H shielding of H shielding” of H 
LiH 39.52 on -14.95 -13.8 +0. 5° 
HF 110.674 110. 9© -87.30 -79.6 +0.5f 





* All results expressed in parts per million. 


Calculated from spin-rotation interaction constants. ! 


d 


Cy=+t8 +1 ke/sec, L. Wharton, L. P. Gold, and W. Klemperer (private communication). 
J. Kolker and M. Karplus obtained 110.68 ppm using more exact values for the molecular integrals. These 


authors are also engaged in calculating the paramagnetic shielding in these molecules with a multiterm variation 


function (private communication). 
f 
(1961). 


©J. N. Pinkerton and C. H. Anderson, Bull. Am. Phys. Soc. 6, 281 (1961). 
Cy=+t71 43 ke/sec, M. R. Baker, H. M. Nelson, J. A. Leavitt, and N. F. Ramsey, Phys. Rev. 121, 807 


representation of C.,, i.e., the same excited states contribute to each term. Hence, we find for the 


magnetic screening constant’ 


LUM Mea Tay lOMOlU2 4/74 "1 


(1) 





1 1 
M4 ger O12 /r 4 10+ 303 


We have evaluated the diamagnetic and para- 
magnetic shielding constants for the proton in 
LiH and HF from Eq. (1).* The results, given 
in Table I, are in reasonable agreement with 
experiment. Applications to the chemical shield- 
ing and diamagnetic susceptibility in other dia- 
tomic systems, including the anomalous proton 
and fluorine shielding in transition metal com- 
plexes, are in progress. Details of these calcu- 
lations will be published later. 

We are indebted to Mr. Russell M. Pitzer and 
Professor M. P. Barnett for much helpful dis- 
cussion and assistance in evaluation of the in- 
tegrals. The calculations were done in part at 
the Massachusetts Institute of Technology Compu- 
tation Center. Acknowledgement is accorded the 
Office of Naval Research for financial aid to one 
of us (C.W.K.). 





'N, F. Ramsey, Phys. Rev. 77, 567 (1950); 78, 699 
(1950); 83, 540 (1951); 86, 243 (1952). 

*For a summary, see J. G. Powles, Reports on 
Progess in Physics (The Physical Society, London, 
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1959), Vol. 22, p. 433. 

%.. Cc. Snyder and R. G. Parr, J. Chem. Phys. 34, 
837 (1961). 

‘4 variational method for eliminating the excited 
state sum in the expression for the diamagnetic sus- 
ceptibility has been proposed by T. K. Rebane, J. 
Exptl. Theoret. Phys. (U.S.S.R.) 38, 963 (1960) 
(translation: Soviet Phys. -JETP 11, 694 (1960)]. 

M. Baker has calculated the diamagnetic susceptibility 
of molecular hydrogen with this procedure and has 
proposed its extension to magnetic shielding (private 
communication). 

5Hartree atomic units (h=1, m=1, e=1) are used 
throughout. 

sy. Schiff, Quantum Mechanics (McGraw-Hill Book 
Company, New York, 1955). For limitations, see 
A. D. McLachlan, J. Chem. Phys. 32, 1263 (1960). 

"K. Ito, J. Am. Chem. Soc. 80, 3502 (1958), has 
also cancelled average excitation energies but in a 
different manner. His method, which was shown later 
[Y. Kurita and K. Ito, J. Am. Chem. Soc. 82, 296 
(1960)] to be equivalent to a variational calculation, 
amounts to placing the origin of the vector potential 
at the center of electronic charge and assuming the 
excited state sum is zero for that origin. 

8Wave functions taken from B. J. Ransil, Revs. 
Modern Phys. 32, 245 (1960). 


















= ass 2h ens 22 


—_. 


rn 


~~ = of rm 6) © = © 89 Cf FF — wr SOOO | a Oo st lle COU, 


Cag ?, a en en) en?) en ee ee 


—=——"-— + O& Ons OO & xX 





961 


(1) 


ter 








VoLUME 7, NUMBER 1 


PHYSICAL REVIEW LETTERS 





Jury 1, 1961 





FIELD-SWEPT MASER OSCILLATION 


James C. Kemp 


Department of Physics, University of Oregon, Eugene, Oregon 
(Received May 1, 1961) 


Recently Singer and Wang’ (SW) have put forth 
an analysis which purportedly describes the be- 
havior of a very wide variety of maser oscilla- 
tors. An experimental verification was offered, 
based on a pulsed maser oscillation waveform 
from the present writer’s work.’ 

The purposes of the present note are as follows: 
(1) To show that our waveforms, on one of which 
SW based their Fig. 1, show large qualitative and 
quantitative discrepancies from their predicted 
radiation envelope structure. For this purpose 
we will first discuss their main Eq. (8) and show 
how its solution (SW Fig. 2) differs from our 
waveforms. Our reason for the discrepancies 
is that our oscillator was of the “field-swept” 
type, in which emission occurs as the Zeeman 
field value H, is swept through resonance, giving 
a radiation envelope determined largely by the H, 
sweep rate and showing a modulation due, we be- 
lieve, to a mechanism unrelated to that invoked 
by SW. (2) To summarize our own analysis? of 
the oscillator behavior, which we feel conforms 
adequately with our waveforms. In particular 
our asymptotic form [Eq. (2) below] for relative- 
ly fast sweep and weak oscillation will be seen to 
agree quantitatively within experimental limits, 
and qualitatively in all details, with an appropriate 
sample (Fig. 3) from a series of oscillograms.? 
We hope this note will help explain waveforms seen 
in other paramagnetic maser oscillators involving 
H, sweep. 

As for the nature of the SW solution, the nu- 
tation-angle substitution’ 6 = f ofahuH (t’)dt’ in 
their Eq. (8), with p£ +uH and N, constant, the 
relevant case, gives 6 +7 16+(r7,) “1 sind =0, the 
simple damped-pendulum equation, where T=2Q/w 
= cavity ringing time and 1,= radiation damping 
time.* This writer noted previously® (see also 
Yariv*) that even with 6(0) ~7 (initial inverted 
state), in the pertinent case the radiation ampli- 
tude H =y@ is given adequately by exp( -t/7) 

x sin|t/(t7,)*] after the first “cycle” or two. The 
resulting emission would have mean duration of 
order t and would show roughly uniform 100% 
modulation with period of order (r7,,)%*/2n. In 
our experiments? 7 = 0.6 usec. Our observed 
durations depended on the sweep rate, and 
ranged? from 7 psec (fastest sweep) to 500 psec 
(slowest sweep), with representative values 500, 
15, and 90 usec, respectively, for Figs. 1, 2, 





FIG, 1. Field-swept maser oscillation, displaying 
emission envelope |H(t)|. Slow sweep, dH)/dt = 800 
gauss/sec. 6=0,001; 7,=0.1lusec. Trace length 
200usec. Fractional Zeeman energy release ~70%. 
Shown here is the rising portion of the envelope, which 
decays to the right. The thickness of the trace in the 
upper part is not noise but unresolved rapid modulation. 
The arrow is believed in the vicinity of passage through 
the line-center resonance (¢ =0), though this point was 
not localized experimentally. Equation (2) of text is not 
applicable. 





FIG, 2. Field-swept emission. dH,/dt=2x10‘ gauss/ 
sec. 6 =0.03; T,=0. lysec. Trace length 30usec. Ow- 
ing to the fast sweep, in spite of about 50% Zeeman en- 
ergy release the modulation pattern is surprisingly like 
that of Eq. (2) except for average shallowness. Period 
of fifth modulation cycle T; = 1. 2usec (theor) = 0. 8 usec 
(exp). 





FIG. 3. Field-swept emission. dH)/dt=4x10* gauss/ 
sec. £20.03; T, =0.6 usec. Trace length 200 usec. 
Energy release ~ 30%. 7T;=2.8usec (theor) =4.0puséc 
(exp). Equation (2) well describes this plus seven other 
oscillograms in reference 2, with moderate to rapid 


sweep and not too strong oscillation. 
21 
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and 3 shown in this note. Our waveforms show, 
rather than uniform 100% modulation, a modu- 
lation with frequency increasing indefinitely with 
time, after sweep through resonance, with a per- 
cent depth that depends strikingly on the sweep 
rate. 

The possibility of observing two-level paramag- 
netic maser oscillation showing the pendulum - 
like behavior® expected from a static-H, anal- 
ysis’’* was investigated by the writer. We con- 
cluded that, when one induces oscillation by the 
common procedure of sweeping the field H, into 
or through resonance, at least in our apparatus 
it was impossible to establish oscillation condi- 
tions “suddenly” enough to allow emission to ac- 
cur with essentially static parameters (the field 
H,). Emission with H, varying in time linearly 
through resonance was then studied systemati- 
cally, and was the subject of our work,” which 
we now summarize. 

Our theoretical treatment? of the field-swept 
case starts with the coupling equation between 
a complex radiation field amplitude H and a com- 
plex transverse magnetization M. For sweep 
through resonance rapid enough to allow only a 
small fractional Zeeman energy release, such 
that M, = -Mo holds, the spin-dynamic equations 
are soluble explicitly for arbitrary H(t) and AH,(t) 
=H,(t)-w/y. We choose AH,=y “at, where a@ is 
a constant sweep rate. We assume a Lorentz 
line of width (T,*)~ and, for convenience,’ that 
T,*=tT. One gets 


aH (2 dH if1 1 
a—aof tituleawte ef =0 
on(? iat) a s(3 ‘. iat) H . (1) 


which is readily transformable to Weber’s equa- 
tion.* The pertinent solution, in terms of a Web- 
er’s function of imaginary order, is 

H=K exp/(iat? /4 -t/1)D; j,l-C +i)(a@/2)"*t], 


with K=const. From asymptotic forms?*® we 
find, for moderate to rapid sweep, 


IH()|=A(|1 + gang; sin(dat? + 6)], (2) 


where ¢ varies slowly and 8 =a7T, is a dimen- 
sionless sweep rate. This form described the ra- 
diation envelope for 8 21/87 =0.04 and t2(27/a)¥?. 
|H(t)| rises monotonically during the interval 
about passage through resonance (t=0). A(t), 
giving the gross shape, rises to simple max- 
imum at roughly ¢, =8~*(27,,/7)“°r, and then de- 
cays, asymptotically like exp(t/r -1/8at?). The 
sine term gives an amplitude modulation, with 
fractional depth about 100% (for 8 ~ 1/87) initially, 
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but decreasing in percentage as 1/t, and increas- 
ing indefinitely in frequency with ¢. 


The modulation mechanism described by (1) and 


(2) inherently involves the time-varying H,, which 
continually tends to throw the magnetization into 
the wrong phase, with respect to the radiation 
field, for oscillation buildup. Modulation due to 
this mechanism should vanish for vanishingly 
slow sweep.°® 

Figures 1-3 are from the writer’s experimental 
field-swept oscillator,” and display |H(t)!. Rele- 
vant measured parameters were: T=0.6usec; 
To*=0.8 usec. tT, and dHg/dt=y ‘a were adjust- 
able, respectively, over the ranges 0.1-0.8 usec 
and 0.8-40 kilogauss/sec, giving a 8 range 0.001- 
0.06. Equation (2) should apply best to Fig. 3. 
The linearization Mz = -Mog of Eq. (1) holds rea- 
sonably for Fig. 3 but marginally for Fig. 2. In 
spite of small 8, and nonlinearity, even Fig. 1 
has some features in common with Eq. (2). Mod- 
ulation appears most marked just when Eq. (2) 
applies. The predicted period T,, =(7/na)” from 
(2), of the mth modulation period to appear, is 
correct within 50% for T, in Figs. 2 and 3. For 
the expected emission duration T, a useful rough 
estimate? is 2(rr,)*B, yielding reasonable values 
500, 20, and 50usec for Figs. 1, 2, and 3, re- 
spectively. 

We have shown that the maser oscillator emis- 
sion envelope had in our case a structure depen- 
dent almost entirely on the field sweep. Though 
a simple pendulum-like energy exchange between 
cavity and emitting sample, described by the 
“general analysis” of SW, is a basic process,’ 
it had nothing to do with our case, with either 
fast or slow H, sweep. 





‘J, R. Singer and S, Wang, Phys. Rev. Letters 6, 351 
(1961). 

3. C. Kemp, thesis, University of California, Berk- 
eley, 1960 (unpublished). The relevant portion is avail- 
able as a University of California Report AFOSR-TN- 
60-509, issue 275, E.R.L., Berkeley, 1960 (unpub- 
lished). 

33. C. Kemp, J. Appl. Phys. 30, 1451 (1959). Though 
this paper covered the rather fundamental “pendulum nu- 
tation” mechanism (discussed earlier, without the non- 
linearity, in reference 4), the possible role of Hy sweep 
in actual experiments was not then appreciated. Adding 
line broadening to the static-parameter (static H»)) emis- 
sion problem was discussed in this reference, in refer- 
ence 5, and further in reference 2 (thesis only), and 
gives a possible lengthening of emission duration and a 
modulation frequency which can be smaller on the aver- 
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age but cannot increase with ¢ over the envelope. 

4n. Bloembergen and R. V. Pound, Phys. Rev. 95, 
8 (1954). ." = YM pMUoV (mks units”), where y =mag- 
netogyric ratio, M)/V = available magnetization per unit 
cavity volume. In SW notation, . =47N,ypQ or 
47NeWwQ. 

5See A. Yariv, J. Appl. Phys. 31, 740 (1960), who 
omitted a minor approximation (reference 2 Report, 
Appendix) leading to the pendulum equation, but gave a 
computer solution justifying our statement that the solu- 
tion is quite like an exponentially decaying sinusoid. 

An unpublished waveform of P, F. Chester, P. E. Wag- 
ner, and J. G. Castle, Jr. , Scientific Paper 6-94439-8- 


P4, Westinghouse Research Laboratories, Pittsburgh, 
Pennsylvania, 1958 (unpublished), p. 10, taken without 
Hy sweep, shows the characteristic 100% almost uniform 
modulation, in my knowledge the only observed such case. 

"It happened that T,*/T=1 in our experiments. The 
mathematics can easily accommodate the more general 
case. 

SE. T. Whittaker and G. N. Watson, Modern Analysis 
(Cambridge University Press, New York, 1935), 4th ed., 
pp. 347-351. 

*Sweep-rate dependence of modulation depth was also 
noted by S. Foner, L. R. Momo, and A. Mayer, Phys. 
Rev. Letters 3, 36 (1959). 








HELICITY OF 4» MESONS; MOTT SCATTERING OF POLARIZED MUONS* 


Marcel Bardon, Paolo Franzini,f and Juliet Lee 


Columbia University, New York, New York 
(Received June 13, 1961) 


Current speculations’ on weak interactions 
raise the possibility that the neutrino associated 
with the » meson may be different from the neu- 
trino of beta decay. Thus the helicity of Vv, May 
be different from that of Vs This experiment 
was designed to observe the helicity of the u~ 
meson (and hence that of the associated Du) 
from the reaction: 


T~>u~+v. 1 
mn . (1) 


The method used was to analyze the polarization 
of the muons by Coulomb scattering of a trans- 
versely polarized u~ beam (Mott scattering) which 
predicts a left-right asymmetry due to spin-orbit 
coupling: 


o(8, ) =f(6)[1 +g(@)P,, cos® |. (2) 


Here P,, is the transverse muon polarization, 
along the y axis, and the z axis is in the direction 
of the incoming muon momentum. Explicit cal- 
culations of the energy-dependent terms, /f(@) and 
g(6), were performed by Rawitscher,? and by 
Franklin and Margolis.* Integration of the equa- 
tion for the experimental arrangement used leads 
to a predicted asymmetry, 


(L -R)/(L +R) =+ 0.09 for P, =F1, (3) 
where L and R represent the number of scattering 
events into the left and right counters with respect 
to the plane defined by the » momentum and spin. 
Our results, corrected for accidentals and 
based on ~ 100 hours of counting (~3 x10 muons), 


were L=515 and R=618. This leads to an asym- 
metry of -0.090+ 0.031 for our Py =0.9, from 
which it is concluded that the helicity of the u~ 
meson in reaction (1) is positive (right-handed) 
in agreement with the V-A theory. Thus the anti- 
neutrino in (1) is also right-handed, and hence no 
evidence is deduced from this result for vp #v,. 
This conclusion is in agreement with those ob- 
tained by measuring the knock-on electrons pro- 
duced by high-energy longitudinally polarized yu 
mesons in magnetized iron.**® 

The agreement of predicted and observed asym- 
metry is interesting in that, to our knowledge, 
there are no pre-existing data on polarized elec- 
tron or muon scattering in this momentum trans- 
fer region ((Aqg),y~100 Mev/c). More explicit 
discussion of this and of o(@) will be published 
later. We only note here that the inelastic scat- 
tering contribution is strongly limited by the low 
energy of the outgoing muon. 

The Mott scattering experiment involved the ex- 
traction of a low-energy 7 beam, and the produc- 
tion from this beam of transversely polarized pu 
mesons by decay in flight. A magnetic channel 
was designed to this end, and its front inserted 
into a recessed window in the cyclotron vacuum 
chamber. The resulting beam of 43-Mev 7 mes- 
ons was further moderated by 3 in. of lithium, 
placed at the focus of a collimating magnetic sys- 
tem so as to produce a parallel (+2°) beam of (28 
+ 2.5)-Mev 7 mesons in the decay region. The re- 
sulting muons selected by decay angle had an av- 
erage transverse polarization of 90% (Py =0.9) 
in the plane of decay, and were uniformly dis- 
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tributed in energy from 9 to 33 Mev. This energy 
band was chosen as a compromise between inten- 
sity of backward scatterings and magnitude of 
asymmetry. The experimental arrangement 

(Fig. 1) makes use of the azimuthal symmetry 

of the pion decays in order to maximize the count- 
ing rate and minimize instrumental asymmetries. 
Lead was used as the scatterer and the scattering 
angles accepted were from 105° to 165°. The 
thickness of the lead was chosen for the greatest 
yield of events consistent with a not too substan- 
tial reduction of the asymmetry by plural scat- 
terings. 

The incoming muon is counted by a thin plastic 
counter (No. 3) and strikes a 2-g/cm? thick lead 
sheet. An event occurs when the muon is scat- 
tered back through No. 3 into either the left or 
the right No. 4 counters which are placed radi- 
ally on either side of No. 3. This setup of count- 
ers and lead is repeated 10 times so as to form 
an annular arrangement about the pion beam as 
axis. The small scattering probability (~ 0.6 
x107° per muon), together with a background 
of 4 -mesonic x rays due to muons stopping in 
the lead, makes it necessary to further identify 
events by detecting the decay electrons from the 
scattered muons which are required to stop in 
the No. 4 counters (3 in. thick). Direct passage 
by a particle from the beam through No. 4 and 
then No. 3 is prevented by 3-in. copper shields 
mounted on the edges of the No. 4 counters to- 
ward the beam. 

Figure 2 shows a simplified block diagram of 
the fully transistorized electronics. The No. 2 
counter was put in anticoincidence with No. 1, 
and ten such 12 signals were obtained from an 
emitter follower fan-out system. Each was then 
put in coincidence with a signal from one of the 
ten No. 3 counters; the 123 then representing a 
~ meson which reached a particular set in the 
wheel. Each of these signals was sent to two 
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FIG. 1. Experimental 
+ _ setup showing disposition 
INTO CENTER of beam counter No. 1, 


SS collimating counter No. 2, 
and wheel arrangement of 
LS ten sets of counters and 





lead scatterers. 


emitter followers, and, as shown in Fig. 2, was 
counted in coincidence with a pulse from an ad- : 
jacent No. 4 counter through a double-pole, 
double-throw relay. This device permitted the 
switching back and forth at regular intervals of 


the whole electronic chain following the 123 co- 
incidence, thus eliminating some possible instru- 
mental asymmetries. The 1234 coincidence out- 
put was used to open a 3-usec gate during which 
a pulse in No. 4, due to an electron from the de- 
cay of a stopped » meson, would count as an 
event. Single counts in No. 4 gave accidental 
events which were about one-tenth of the true 
events. This low accidental rate was made pos- 
sible by the use of the Nevis vibrating target,° 
which gave an average duty cycle of 3. 
Instrumental asymmetries were rendered neg- 
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ligible by the following steps: 

(i) During the 95 hours of counting, the elec- 
tronic channels for left and right scattering were 
switched every hour. 

(ii) The arrangement is such that each of the 
No. 4 counters is a “left” counter for one set 
and a “right” counter for the adjacent set. 

(iii) The counters were aligned on the wheel 
to within +0.016 in. around the 3-ft diameter. 
The repetitive pattern also serves to average 
out any residual asymmetry. 

(iv) The wheel assembly was rotated about 
the 7 beam by one-fifth turn (incoherently) every 
two hours to smooth the small deviations from 
axial symmetry (beam shape and aim). 

(v) Finally, a 7* beam was scattered into the 
counters by 6 g/cm? of lead placed at No. 2 col- 
limator. The 7* mesons simulated true events 
by scattering back from the lead, stopping in 
No. 4 and giving a delayed electron via the 1* + 
u++e* chain. This yielded an inherent asym- 
metry (for almost isotropic scattering from Pb) 
of +0.008 + 0.019. 

The authors wish to express their gratitude to 
Professor Leon M. Lederman for his continuous 





guidance and encouragement; and also to thank 
Uriel Nauenberg, Michael Tannenbaum, and the 
staff of the Nevis Cyclotron Laboratory for their 
assistance at various stages of the preparation 
of the experiment. 
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RELATIVE INELASTICITY AND ANISOTROPY OF PROTON INTERACTIONS AT 9 AND 23.5 Gev 


E. M. Friedlander, M. Marcu, and M. Spirchez 


Cosmic-Ray Laboratory, Institute of Atomic Physics, Bucharest, Rumania 
(Received June 15, 1961) 


Recent observations of cosmic-ray jets’ have 
brought evidence for a decrease of the inelastic- 
ity K of nucleon-nucleon collisions with increas- 
ing energy of the primary, together with an in- 
crease of the forward-backward peaking of the 
c.m. angular distribution of the secondaries. 

In view of the obvious uncertainties connected 
with the estimation of primary energy and/or 
energy and identity of the created particles, we 
thought it desirable to check these trends under 
as far as possible controlled conditions at the 
highest available accelerator energies. We have 
measured K by two independent methods and 
found it to decrease significantly in the energy 
interval 9-25 Gev. 

Electron pairs from 7°-2y decay were detected 
by upstream scanning? along relativistic tracks: 
(1) in a stack of NIKFI-R pellicles, used in an 


earlier investigation,® irradiated in the internal 
beam of the Dubna proton-synchrotron circulating 
at 9-Gev total energy, and (2) in two Ilford G-5 
pellicles irradiated in the external (scattered) 
momentum-analyzed beam of the CERN proton 
synchrotron (average momentum ~ 23.5 Gev/c). 
Only tracks with a projected length per plate > 2 
mm and a projected angle < 30° with the beam 
direction were accepted. For each detected pair 
the projected opening angle a and the angle 6 of 
the pair bisector with the beam direction were 
recorded. Since at least one electron track in 
each pair was essentially flat, the quantity 

(a) is proportional to the average 7° energy 
within the accepted solid angle interval. In order 
to obviate the well-known difficulties connected 
with the estimation of the proportionality factor,* 
we preferred to compare directly (a) from the 
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two stacks, since 


-1 
Ko3.5_ 9 6%) 93.5 ("593.5 


—— £ : 
Ky 23.5 (a dg (ng 





(1) 


the inequality sign in (1) being due to the stronger 


collimation of the secondaries at the higher en- 
ergy. Using the experimental results on (a~!) 
given in Table I and the charged-pion multiplic- 
ities detected in along-the-track scans,”’® we 
obtain 


Kys.5/Ky < 0.70 + 0.08, (2) 


The quantities (@~'sin@), which are a meas- 
ure of the average transverse momentum of the 
neutral pions, are also given in Table I. The ab- 
sence of any significant difference in the trans- 
verse momenta, as well as the practical inde- 
pendence of the transverse momentum on emis- 
sion angle, allows an independent check of the 


above result by means of the charged secondaries. 


Meson showers of at least three relativistic 
prongs and N; < 4 were collected by an area scan, 
and the angular distribution of the shower tracks 
was measured. The condition N, <4 selected 
light target nuclei® and provided us with a sample 
enriched in proton-nucleon collisions. No effort 
was made to detect Np =0 stars in this run. 

Provided py is independent of emission angle, 
we have 

Ko3.5 9 (€8¢8.)99 54103 5 


K, 23.5 (esc), (n').” (3) 





9 


where @, are the emission angles of the relativ- 
istic secondaries with respect to the incoming 
proton, and (n,’) is the average multiplicity of 
shower tracks in the analyzed showers (7, > 3). 


Table I, Comparative data from measurements at 9 
and 23.5 Gev. 








Stack 9 Gev 23.5 Gev 
No. of electron pairs 504 102 
(a~!) 251411  273+27 
(a-! sind) 37 +4 35 +4 
No. of stars with n >3, N, <4 616 67 
No. of tracks 2281 441 


a’) 4.33+0.05 4.92+0.23 
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FIG, 1. Integral normalized distributions of cscé; 
for 9-Gev meson showers (squared points) and 23,.5- 
Gev meson showers (circles) with n, #3, N,<4. 


The (normalized) integral distributions of cscé, 
are shown in Fig. 1. The pertinent numerical 
data are also given in Table I and yield 


Ky5,5/K, = 0.67 + 0.06, (4) 


a result quite consistent with that from the 7° 
measurement. Using the upper limit for the 
inelasticity at 9 Gev, deduced in reference 4, 
we obtain 


K,5.5 < 0.23 + 0.06. (5) 


In order to check the peaking of the c.m. angu- 
lar distribution, an anisotropy parameter, 


n = ($7)¥( olin cote | -(1n coté -) [mm -1))™, 


(6) 


was computed for every shower with no relativis- 
tic prong in the backward hemisphere. It can be 
shown that if the quantities In cot@, belong to a 
normal population of variance o”, then 


{n) =0, (7) 
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—s and : played by peripheral collisions, involving “almost 
o =0.75nn -¥ (8) real” pions in the meson clouds of the collision 
ins a partners.°~ The implications of these results 
will be discussed quantitatively in a forthcoming 
| The distributions of the 7 values for 9-Gev and paper, based upon the model discussed in refer - 
23.5-Gev showers are shown in Fig. 2 normalized ences 10 and 11. 
to the same area. The marked shift towards high- We are greatly indebted to Professor V. F. 
er 7 values at the higher energy is proved to be Weisskopf and Dr. W. O. Lock for making avail- 
significant by a Pearson test (x? =20.5 with 6 de- able the plates exposed at the CERN proton syn- 
grees of freedom). The average values are chrotron, to Professor J. Auslander for his sup- 
port in this matter and for useful discussions, to H. 
(1) 23.5= 0.54 + 0.02, Totea for his assistance in processing the angular 
\ data, and to our scanner team for their whole- 
i (n),=0.42 + 0.01, (9) hearted cooperation. 
and prove that both distributions are significantly 
peaked (0.39 would be expected for isotropy). ‘See, e.g., D. H. Perkins, Progress in Cosmic-Ray 
The above results, which are consistent with Physics, edited by J. G. Wilson (Interscience Publish- 
a ‘ ers, New York, 1960), Vol. 5; also L. Montanet, J. A. 
recent aren evideace from counter experi- Newth, G. Petrucci, R. A. Salmeron, and A. Zichichi, 
ments,‘** are suggestive of an increasing part Nuovo cimento 17, 166 (1960); S. Lal, Y. Pal, and A. 
Raghavan (to be published). 
*p. I. King, Phys. Rev. 109, 1344 (1958). 
aN 3E. M. Friedlander, Nuovo cimento 14, 796 (1959). 
ae an! 4G. L. Bayatyan, I. M. Gramenitzkii, A. A. Nomofi- 
2scé 16 — lov, M. S. Podgoretzkii, and E. Skrzypezak, Zhur. 
\ ' Eksp. i Teoret. Fiz. 36, 690 (1959) [translation: So- 
we 2 4 asa, Ot viet Phys. —JETP 36(9), 483 (1959)]. 
s ' : 25 Bev I, M. Gramenitzkil, M. Ia. Danysh, V. B. Liubimov, 
, mo M. I. Podgoretzkil, and D. Tuvdendorzh, Zhur. Eksp. i 
' | Teoret. Fiz. 35, 552 (1958) (translation: Soviet Phys. — 
| L JETP 35(8), 381 (1959)]. 
R . 6G. Cvijanovich, B. Dayton, P. Egli, B. Klaiber, 
Os 4 ' W. Koch, M. Nikolic, R. Schneeberger, H. Winzeler, 
| Sy J. C. Combe, W. M. Gibson, W. D. Lock, M. Schnee- 
: L berger, and G. Vanderhaeghe (to be published). 
os ; 'G. Cocconi, A. N. Diddens, E. Lillethun, and A. M. 
: Wetherell, Phys. Rev. Letters 6, 231 (1961). 
(4) 4 ' "a 8M. Fidecaro, G. Gatti, G. Giacomelli, W. A. Love, 
: ; W. C. Middelkoop, and T. Yamagata, Nuovo cimento 
2r — 11, 382 (1961). 
| ae es his. coger ~ 9N, Yajima, S. Takagi, and K. Kobayakawa, Progr. 
: 0 0.2 0.4 0.6 0.8 1.0 12 . Theenet. Phys. (Kyoto) 24, 59 (1960). 
E. M. Friedlander, Phys. Rev. Letters 5, 212 
FIG, 2. Differential distributions of the anisotropy (1960). 
} parameter 7 in the same showers. Dashed line at 9 Gev, ‘lf, Salzmann and G. Salzmann, Phys. Rev. Letters 5, 
(5) H full line at 23.5 Gev. The two histograms are normal- 377 (1961). 


ized to the same area. 125. Dremin and D. Cernavsky (to be published). 
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SCATTERING OF 915-Mev a PARTICLES FROM CARBON AND HELIUM: 
DIRECT EVIDENCE FOR a-PARTICLE CLUSTERING IN NUCLEI" 





Terence J. Gooding and George Igo 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received June 2, 1961) 


We have obtained direct evidence for a -particle 
clustering in C™ from the (a, 2a) reaction at 915 
Mev. Measurements of the energies of the two 
outgoing a particles and the angular correlation 
between them show that the reaction proceeds via 
a direct collision between the incident a particle 
and an a-particle cluster in the nucleus. An esti- 
mate of the preformation probability has been ob- 
tained by comparing the carbon and helium dif- 
ferential cross sections. 

Figure 1(a) is a diagram of the experimental 
arrangement. The 915-Mev a-particle beam 
from the Berkeley synchrocyclotron was focused 
at the center of an 18-in. diam scattering chamber 
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FIG. 1. (a) Diagram of the experimental arrange- 
ment. (b) Block diagram of the electronics. 
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to a spot 1 in. wide by 3/4 in. high. The beam 
intensity was about 3x10® @ particles per sec and 
was monitored by an 8-in. diam ionization chamber 
The duty cycle was increased from about 60 psec 
per burst to 12 msec per burst by means of an_ | 
extra peripheral dee. Measurements were made | 
first with the chamber filled to 65 psi of helium | 
and then with a thin graphite target (~ 260 mg/cm’ | 
for 6; =17 deg and ~600 mg/cm? for 67 = 26 deg) 
placed at the center of the scattering chamber. 

Particles scattered from the target were de- 
tected by counters placed at angles 6, and 6p 
with respect to the a-particle beam. The parti- 
cles were identified by measuring their total en- 
ergy, E, and their rate of energy loss, dE/dx. 
The dE/dX counters were disks of plastic scin- 
tillator, 1.0 in. diameter by 0.1 in. thick, and the 
E counters were cylinders of plastic scintillator, 
3.5 in. diameter by 11 in. long. The response of 
the counters to a particles and protons of various 
energies was determined by placing the counters, 
preceded by an appropriate absorber, in the di- 
rect a-particle and proton beams. The response 
was in reasonable agreement with calculations 
based on the response of plastic scintillator to 
lower energy protons, deuterons, and electrons.’ 
The nuclear absorption in the E counters was 
about 35% and was determined in the direct beam 
by measuring the low-energy tail of the pulse- 
height spectrum. 

A block diagram of the electronics is shown in 
Fig. 1(b). The fast-coincidence resolving time 
was 5 nsec and was limited by the flight-time dif- 
ference between high- and low-energy a particles 
leaving the target. The particle-selection system 
consisted of complex arrangements of discrimina- 
tors and coincidence circuits acting on linear ol 
puts from the dE/dX and E counters. In this way, | 
we were able to sort doubly charged particles j 
(He® and He*) from singly charged particles with | 
a high degree of reliability. Outputs from the 
particle selection systems were fed into a slow- 
coincidence unit which was used to start a mag- 
netic recorder. With the helium target, the fast- 
coincidence counting rate was typically one per 
second, whereas the (a, 2a) coincidence rate was 
only one per minute. Outputs from the E counters 
were pulse-height analyzed in the magnetic re- 
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corder. In addition, information concerning the 
energy range in which each a particle fell was 
recorded. The magnetic tape was subsequently 
analyzed using an IBM-7090 computer. 

The differential cross section d*0/dQ,dQ, for 
scattering of the a-particle beam by helium and 
carbon was measured at 6; =17, 26, and 45 deg, 
with 67 + @p= 87 deg (the separation angle for 
elastic scattering at 915 Mev). At each angle, 
we measured the energy spectrum E, and the 
spectrum of the sum of the energies E; +Ep. 

The spectra from helium at 17 deg are shown 
in Fig. 2. The peak at 815 Mev in E, and 915 
Mev in (E,; +E R) corresponds to the elastically 
scattered group. The second peak at E = 650 
Mev and (E; +Ep)= 750 Mev is thought to be due 
to stripping of the forward-going a particle into 
a He® plus a neutron. The relative intensities of 
these two peaks are misleading, since the elastic 
peak is sharply correlated at these angles, where- 
as the inelastic group would be expected to have 
a broader distribution. The third peak may be 
due to an energy cutoff in the experiment at around 
350 Mev. The energy calibration throughout these 
measurements is uncertain to +50 Mev. At 26 deg 
the spectrum looked essentially the same but the 
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FIG. 2. Scattering from helium at 67 =17 deg, 


67 +8p=87 deg. (a) Summed energy spectrum E; + Ep; 
(b) energy spectrum Ey. The energy calibration is 
uncertain to +50 Mev. 


cross section was reduced. The cross section at 
45 deg was too low to be measured. 

The carbon data at 26 deg are shown in Fig. 3. 
The summed spectrum again shows two peaks, 
but the relative intensity of the elastically scat- 
tered group appears reduced. This is because 
the angular correlation is broadened by the mo- 
mentum distribution of the a particles in the car- 
bon nucleus. The peak at 750 Mev in the energy 
spectrum E, for events with E; +Ep > 850 (i.e., 
the elastically scattered group) shows clearly that 
the reaction is proceeding via a direct collision 
of the a particle with an a cluster in the nucleus. 
If, for instance, the initial collision was between 
the incident a particle and a deuteron, then the 
energy of the a particle scattered at 26 deg would 
be about 500 Mev. An a-particle—nucleon col- 
lision producing an a particle at 26 deg is kine- 
matically forbidden. 

Further evidence in support of this conclusion 
is given in Fig. 4 which shows the summed spec- 
tra for 6; =17 deg and 67 + @p equal to (a) 87 deg 
and (b) 57 deg. The intensity of the quasi-elas- 
tically scattered group drops by about a factor of 
four when the separation angle is moved 30 deg 
from the kinematic angle of 87 deg. Consequently, 
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FIG. 3. Scattering from carbon at 67 =26 deg, 
6; + Op =87 deg. (a) Summed energy spectrum E,; + Ep; 
(b) energy spectrum E;, where E, +Ep >850 Mev. 
The energy calibration is uncertain to +50 Mev. 
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FIG, 4. The carbon summed energy spectrum at 
6, =17 deg and (a) 07 +p =87 deg; (b) 6; +9p=57 deg. 
The energy calibration is uncertain to +50 Mev. 


the recoil a particle is correlated to within (6+ 3)% 


of the total available solid angle. 
In order to obtain the probability of a -particle 
clustering in carbon, the differential cross sec- 


tions for elastic scattering from helium and quasi- 


elastic scattering from carbon have been calcula- 
ted at 17, 26, and 45 deg and are shown in Table I. 
The differential cross section for helium has been 
corrected for the finite size of the beam, the fi- 
nite extent of the helium gas target, and nuclear 
absorption in the counters. To obtain the cross 
section for carbon, one must integrate over the 
solid angle into which the correlated a particle 
can be scattered. 

To perform this integration, we have assumed 
that the angular correlation follows a triangular 








Table I. The differential cross sections for elastic 
scattering from helium and quasi-elastic scattering 
from carbon. 











“ da/dQ (ub/sr) da/dQ (carbon) 
(deg) Helium Carbon do/dQ (helium) 

17 165 +35 287 +150 1.7 26.9 

26 22 +11 86 +50 3.9 +2.9 

45 =“ =e 








distribution through the two experimental points, 
that the angular correlation at 6; = 26 deg is the 
same as at 6, =17 deg, and that the width of the 
distribution out of the plane of scattering is the 
same as that in the plane of scattering. The large 
errors on the cross sections quoted in Table I 
are primarily due to the inadequacy of the data on 
the angular -correlation function. 

The ratio of the (a, 2a) differential cross section 
for carbon and helium is 1.74+0.9 at 17 deg and 
3.9+ 2.9 at 26 deg. According to the extreme a- 
particle model for carbon, we would expect this 
ratio to be about 3. 

Further measurements are in progress on the 
angular -correlation function and the differential 
cross section for the (a,2a) reaction for other 
nuclei. 

We wish to thank Charles A. Corum for design- 
ing the scattering chamber, Melvin Brown and 
Ellwood S. Douglas for assistance with the elec- 
tronics, Dr. Paul Concus for the computer pro- 
gram, and the crew of the 184-in. cyclotron for 
their cooperation. 





*This work was supported by the U. S. Atomic 
Energy Commission. 

'T, J. Gooding and H. G. Pugh, Nuclear Instr. and 
Methods 7, 189 (1960). 
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NUCLEAR DE-EXCITATION FOLLOWING MUON CAPTURE 
AND THE BOUND MUON DECAY ANOMALY* 


Frank Chilton 


University of Washington, Seattle, Washington 
(Received June 9, 1961) 


Experiments have been performed** to meas- 
ure the rate of negative muon decay, 


+e +v4+D, 
when the muon is bound in the Coulomb field of 
| anucleus. Figure 1 presents the data of Yovan- 
| ovitch’ for the ratio of the bound decay rate to 


| the free decay rate (rate for u* decay or rate 
for muons bound to a low-Z nucleus), 


- + 
R=A, /A,- 


The curve is a theoretical curve.* While the 
theoretical predictions are always less than 
R=1 the experimental results rise well above 
R=1 in the iron region, yet return quickly close 
to theory for ,,Zn, and fall well below the theo- 
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FIG, 1. Data of Yovanovitch (reference 1) on bound 
muon decay. 





retical curve for very large Z. The presence 
of this anomaly has been verified in detail by 
combining measurements by Lederman and 
Weinrich? of Aq /M¢ and measurements of A; 5 
the total disappearance rate. These points are 
not included in Fig. 1 simply to avoid clutter. 

The other mode of disappearance of negative 
muons in matter is muon capture for which the 
basic reaction is 


LL +pon+v. 


To determine whether the bound muon decay 
anomaly is really due to detection of y rays 
associated with muon capture rather than elec- 
trons from muon decay, we can take the differ- 
ence between the experimental and theoretical 
values for R and divide by A,/Ag*. Table I 
presents the values of (Rexp -Rtp)Ag’ /A¢ calcu- 
lated using experimentally determined capture 
rates.° The near constancy of the values in 
Table I makes it appear that the anomaly is due 
to the detection of 1% of the muon capture events 
and that this sensitivity to muon capture myster- 
iously stops at Zn. The initially sharp Z depend- 
ence of the anomaly is just the Z dependence of 
muon capture rates ~Z*. 

To understand how background due to muon 
capture could appear in such careful experiments, 
let us consider some of the nuclear processes 


‘Table I. Proportionality of the bound muon decay 
anomaly to the capture rate. 





+ 
(R apa d /n., 





Z Element («« 10?) 
20 Ca 0.4 40.5 
22 Ti 0.7 +0.5 
23 Vv 1.5 +0.4 
25 Mn 1.0 +0.4 
26 Fe 1.4 +0.3 
26 Fe 1.7 +2. 0% 
27 Co 1.3 20.3 
28 Ni 0.8 40.3 
29 Cu 1.2 +0.9% 
30 Zn -0.1 40.2 














“From measurements of Aq /Az by Lederman and 
Weinrich (reference 2). 
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accompanying muon capture in complex nuclei. 
For capture by a free proton the recoil neutron 
has an energy of 5 Mev. In complex nuclei the 
average nuclear excitation is 10 to 15 Mev de- 
pending on the nuclear model.** All but about 

2 Mev is imparted to the neutron. 

When the neutron excitation is greater than the 
neutron separation energy, neutron emission is 
usually most likely. From the experimental cross 
sections for (x,y) reactions, which are typically 
a few millibarns,°® and from experimental cross 
sections for neutron scattering, which are typi- 
cally a few barns, the ratio of the photon to neu- 
tron emission rates is 


r /T ~1073, 
Y n 


when the excitation is a few Mev above the neutron 
separation energy. In particular this indicates 
that in this energy region the emission of photons 
with £,,> 10 Mev is unlikely.*°™ Other properties 
of neutron emission have been considered in con- 
nection with the problem of the neutron asymme- 
try from the capture of polarized muons.”*”*!° 
When the neutron excitation energy is less than 
the separation energy, only de-excitation by pho- 
ton emission is possible. Note that this situation 
can occur as the second stage of compound nu- 
clear decay even when the initial excitation is 
greater than the separation energy. Also note 
that compound nuclear processes are expected 
to prevail over direct neutron emission.** While 
it would certainly be difficult to calculate the 
photon spectrum that results from muon capture, 
a great deal of qualitative information can be ob- 
tained from experimental data on the y-ray spec- 
tra resulting from thermal neutron capture.” 
This corresponds to an excitation close to the 
separation energy, the highest excitation result- 





ing only in photon emission. We have already 
seen that this situation occurs with sufficient 
probability. 

The first thing to be learned from thermal 
neutron capture is an estimate of the mean mul- 
tiplicity of photons. This varies between 2.5 to 
4 photons per capture for the medium-weight 
nuclei, depending on the nucleus involved. The 
mean multiplicity for muon capture might be ex- 
pected to be slightly higher for each nucleus than 
its thermal neutron counterpart since the excita- 
tion takes place in such a different way and be- 
cause the more favored states in muon capture 
differ from those of photon emission. The “al- 
lowed” and “first forbidden” transitions in muon 
capture correspond to 0+, 1+ and 0-, 1-, 2- tran- 
sitions, respectively.** However, since k,R~2.5 
for medium nuclei, a change of two units of for- 
biddenness should make a difference of no more 
than about 10~ in the transition probability.”® 
The most important consideration is which states 
are available. 

Another thing to be learned from neutron cap- 
ture is the presence of strong direct transitions 
in certain regions of the periodic table. In the 
iron region there are strong direct s + 2p,, tran- 
sitions. An example, the thermal neutron capture 
spectrum of Fe’, is shown in Fig. 2. The E,, = 7.85 
Mev, 40% frequency peak corresponds to a tran- 
sition to the 2p,, single-particle state. The small- 
er E, =6.0 Mev peaks correspond to the 2),,. 
state. For reasons discussed below, the bound 
muon decay anomaly seems to be due to y rays 
with 5 Mev < Ey <10 Mev. The experimental 
anomaly in muon decay found in the iron region 
suggests that these same p states seem to play 
an important role in the y-ray transitions that 
cause the anomaly. This is illustrated by the 
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FIG. 2. Thermal neutron 
capture y-ray spectrum of 
Fe (from reference 14). 
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cases 


agNi™?(u -, v)gg7CO™, 
and 
go 20470698 (4, “ ne, 


for which the final nuclei have the p,, states es- 
sentially full.*® Table I indicates that the bound 
muon anomaly is diminished for ,,Ni and gone for 
gon, but not for ,,Co or ,,Cu for which the daugh- 
ter nuclei have empty p,, states. The single- 
particle states which would seem most favorable 
to production of such y rays in muon capture are 
listed in Table II for a few nuclei. This is based 
on the fact that the outermost proton shells are 
the most important in muon capture,’” and on 
considerations of the relative contributions of the 
various states involved.’”*® Configuration mixing 
is apparently responsible for making the p states 
important, which is consistent with evidence’® for 
the applicability of the Nilsson model” in the iron 
region. It is clear from Table II that the yield of 
higher energy y rays resulting from muon capture 
in Ni and Zn will be less than that in the remain- 
der of the iron region. It is also understandable 
why the anomaly disappears at ,,Zn; the 2p,, pro- 
ton shell is being filled but the 2p,, neutron shell 
is already full. However, note that the anomaly 
might still be expected to be present for ,,Cu al- 
though in possibly reduced intensity because the 
shells are partially filled. The data of Lederman 
and Weinrich? are consistent with this (see Table 
I). 

The experimental arrangement of Yovanovitch! 
included a thick target and an electron telescope 
consisting of three thin plastic scintillators with 
two aluminum absorbers sandwiched between 
the scintillators. The absorber thickness was 
equal to the range of a 4-Mev electron. The 
arrangement of Lederman and Weinrich was 
similar.? A 5- to 10-Mev y ray can produce 


either a pair or a Compton electron with en- 
ergy >4 Mev in either the target or the first 
absorber which will be detected by two counters. 
In coincidence, a Compton electron from a low- 
energy y ray can be detected in the remaining 
counter. From the known cross sections” the 
total efficiency for such a high-low combination 
as that above is probably close to 1%, primarily 
because there are several different ways in 
which this event can take place. The discrepancy 
for large Z is evidently due to the different ef- 
fect of thick targets and absorbers on the altered 
electron spectrum for bound muon decay as op- 
posed to the spectrum for free decay.” 

Whether or not the above mechanism is re- 
sponsible for the bound muon decay anomaly 
can be tested experimentally. Use of an exper- 
imental method which absolutely guarantees ex- 
clusion of electrons and photons with E<10 Mev 
would demonstrate whether the anomaly must be 
due to other effects than muon decay. This seems 
to be the case.”* Anderson et al." have measured 
the y-ray spectra from muon capture in Fe and 
Cu. The part of the spectrum with E,,>5 Mev 
does not seem inconsistent with the above ideas. 
The contrast in the y-ray spectra in the iron re- 
gion with that in and after Zn should be particu- 
larly interesting. Whether the high-low combi- 
nation of y rays is likely following muon capture 
could be determined by a coincidence experiment. 
Performing an experiment similar to that of Yo- 
vanovitch’ and Lederman and Weinrich? with an 
electron telescope consisting of many counters 
would be particularly enlightening if all events 
were stored in which, for instance, the first and 
last counters register. Such an experiment 
would permit measurement of bound decay rates 
and examination of the detailed effects of photons 
from muon capture simultaneously. 

It might also be expected that strong direct y- 
ray transitions following muon capture show up 


Table Il, Expected dominant single-particle transitions resulting in y rays with Ey>5 Mev. 








Initial Excited Configuration Final 
Nuclei proton state neutron state mixing neutron state Comment 
- 1 f lg 2 dy 2Py2 
Fes , M 56* { v2 v2 
oFe""(H™ , v)_5Mn 1 fys lew ot 1 fye 
Co™ py» full 
58,60, ,,.— 58, 60* »Py2 
ogNir’ (Me™ , v)gxCo lfy, lgy2 1 fy { Co (py) 
39 Zr 968 ( .— | vag Cue * 2Py2 2dy2 2Py2 Py full 
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in other regions of the periodic table. This 
should be most likely in the light nuclei and near 
closed shells as is the case in thermal neutron 
capture. 

The author would like to thank Professor E. M. 
Henley for many helpful discussions, and Pro- 
fessor H. L. Anderson, Professor J. W. Keuffel, 
and Dr. R. W. Huff for notification of their results 
prior to publication. 
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ANOMALY IN MESON PRODUCTION IN p+d COLLISIONS* 


Norman E. Booth, Alexander Abashian, and Kenneth M. Crowe 


Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received June 12, 1961) 


In an earlier Letter’ we reported the results 
of some preliminary measurements of the mo- 
mentum spectra of He® and H® nuclei produced 
in collisions of high-energy (624- to 743 -Mev) 
protons with deuterium. In addition to the single- 
pion reactions 





He*+ 7° (1a) 
pea] H° +n", (2a) 

we looked for the reactions 
He*+ w® (1b) 
p+d-+| He’+at+77 (1c) 
He® + 7°+ 7°, (1d) 

and 

H?+ wt (2b) 
p+d- H3+ 17+ 7°, (2c) 


where w may be a particle of mass between 1 and 
2.8 pion masses. For reactions resulting in a 








FIG. 1. Experimental 


He’, the two pions (or particle) can be in isotopic 
spin states 0 or 1; if a H® results, only J=1 is 
allowed. We found an anomalous peak in the He® 
spectra which appeared to behave kinematically 
like a particle or resonance of mass approximate - 
ly 310 Mev. At that time we were unable to give 
a definite isotopic spin assignment to the anomaly, 
and we considered a P-wave 7-7 resonance as a 
possible explanation.? 
_ We have since repeated the experiment with a 
new arrangement which enabled us to measure 
both the He® and H® spectra with improved reso- 
lution and accuracy. With the new data we have 
been able to assign an isotopic spin J=0 to the 
anomaly, and subsequently to rule out the P-wave 
m-7 resonance hypothesis. It is possible to ex- 
plain the anomaly by means of a strong S-wave 
m-7 attraction in the J=0 state, with a scattering 
length between 2 and 3 pion Compton wavelengths. 
Figure 1 is a schematic drawing of the experi- 
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mental arrangement. The proton beam extracted 
from the 184-inch cyclotron was passed through 
a gaseous deuterium target operated at liquid 
nitrogen temperature and at about 300 psi. Par- 
ticles produced at 11.8 deg were collimated by a 
system of slits, and focused at infinity by the 
quadrupole @,. Momentum analysis was accom- 
plished by the bending magnets M, and M,. Quad- 
rupole Qe focused the particles at a grid consist- 
ing of six 4-inch-wide counters, each of which 
defined a momentum bite Ap/p of 0.45%. He® and 
H® were selected from other particles by time of 
flight, range, anddE/dx. Backgrounds were meas- 
ured by using hydrogen gas in the target. 

The results of the measurements® of the He* 
and H® spectra at the full proton energy of 743 
Mev are shown in Figs. 2(a) and 2(b). We were 
unable to observe He* of momenta less than 1000 
Mev/c because of their low range. The peaks at 
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FIG, 2. (a) Momentum spectrum of He’® at 11.8 deg, 
laboratory system. (b) Momentum spectrum of H®, 
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1530 Mev/c correspond to Reactions (1a) and (2a) 
at 50 deg c.m. for the heavy particle. The peak 
at 820 Mev/c in the H® spectrum corresponds to 
156 deg c.m. in Reaction (2a). Double-pion pro- 
duction for this laboratory-system angle is kine- 
matically possible between the limits of 910 and 
1440 Mev/c. 

If we assume the anomaly is not due to a pion- 
nucleon interaction, and that charge independence 
holds, we can analyze Reactions (1b), (1c), (1d), 
(2b), and (2c) in terms of the isotopic spin of the 
two pions (or particle). Since the p+d system has 

=4 and J, =4, we can form with He® or H® the 
isotopic spin eigenfunctions having J=4, I, =}4: 


II,°He*, (3a) 
(3)711,2 HS -(3)711,° He’, (3b) 
where, for two pions, 
1,*= @)2(rtn? - 7°n*), (4a) 
1,°= (@)4(x*1~ -1-n*), (4) 
and 
1,°= G)Y*(ata + aa - 20), ee) 


Thus, with H® we have only J =1 production; but 

with He® we have both J=0 and J=1, the amount 

of J=1 being $ as large as in the H® case. To get 

the J=0 spectrum shown in Fig. 3, we drew a 

smooth curve through the H® spectrum, divided 

it by two, and subtracted it from the He® spectrum. 
Before looking carefully into possible explana- 
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FIG, 3. I=0 part of He* spectrum. The dashed 
curve is the phase-space volume fitted to the points 
below 1300 Mev/c. The solid curve is the phase-space 
volume multiplied by the pion-pion enhancement factor 
for a scattering length a,y=2.8h/uc. The experimental 
resolution has been folded into both curves. The wscale 
gives the total energy in the two-pion barycentric system. 
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tions for the bump that appears in the J=0 spec- 
trum, one must know the momentum resolution 
of the experiment. This resolution was calculated 
by taking into account such effects as finite angular 
definition, beam dimensions, image and grid sizes, 
multiple scattering, angular divergences, and 
energy spread of the proton beam. By making 
reasonable assumptions about the energy spread 
of the proton beam, we were able to reproduce 
the shapes and widths of the three single -pion 
production peaks. Figure 4 shows the resolution 
function for He? momenta near 1400 Mev/c. The 
full width at half maximum is 35 Mev/c, with an 
estimated uncertainty of t$ Mev/c. 

As a starting point for comparing the data with 
theory, we computed the Lorentz-invariant phase - 
space volume element in the laboratory system: 


__@p ptf, st ig 
S dpdQ, ws w F 


where 


dp. dp D _ - i 
px APs APs OPs 55 + Pat Ps -P,)5(wgt wat ws - Wz), 
Ws Wa Ws 
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FIG, 4. Effect of subtracting the phase-space volume 


from the J=0, He® data. The solid curve is the computed 
experimental resolution at 1400 Mev/c. 


P, is the momentum of the incoming proton, W, 
is the total energy in the laboratory system, p, 
and w, are the respective momentum and total 
energy of the He*, and w is the total energy in the 
barycentric system of particles 4 and 5, the two 
pions of mass yp. We have assumed the transition 
matrix element to be a constant, and restrict 
ourselves to the relativistically invariant form of 
the volume in phase space. The calculations were 
done for charged and neutral pions and combined 
according to Eq. (4c). With the resolution folded 
in and normalized to the experimental data at mo- 
menta below 1350 Mev/c, ¢g is shown as the 
dashed curve in Fig. 3. 

A promising explanation at the moment is that 
the anomaly is due to a strong S-wave 7-7 inter- 
action* that can be characterized by a scattering 
length. The conditions for the validity of the 
theory of final-state interactions® are: that the 
mechanism of the primary reaction be a short- 
range interaction, that the final-state interaction 
be strong and attractive, and that we consider 
only low relative energies of the two pions. Under 
these conditions, which we will see can be satis- 
fied here, the volume element in phase space 
[Eq. (5)] for a given value of the pion-pion energy 
w is enhanced by a factor proportional to the pion- 
pion scattering cross section at the energy w. To 
obtain the energy dependence of the pion-pion 
cross section we go to Eq. (V.22) of Chew and 
Mandelstam,® where we define a scattering length 
asQ as the amplitude at zero energy. We have 


ve ve 
(fs) coté,° reg In [+ (v+1)], 
4 #50 


+1 1 
(6) 


where 6,° is the S-wave 1-7 phase shift in the J=0 
state and vy is the square of the momentum (in 
pion mass units) in the two-pion barycentric sys- 
tem. This gives an enhancement factor (normal- 
ized to unity at y=0), 


24.6 i on . 
h. 7 (—) info (+1) 


«1 


" 459 () | ? 
(7) 


This is exact for the S-dominant solutions of the 
1-7 equations’ where the 7-7 coupling constant is 
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For the P-dominant solutions® Eq. (7) is still a 
very good approximation, but to evaluate \ we 
must know vp andT’, the position and width param - 
eters, respectively, of the P-wave 7-7 resonance.® 
We obtained good fits to our J/=0 data at momenta 
above 1350 Mev/c with agg between 2.2 and 3.0h/ uc 
with a best-fit value of 2.5%/uc. The solid curve 
of Fig. 3 shows the fit obtained for agg=2.8h/ uc. 
We do not expect the computed curve to fit the ex- 
perimental points below 1350 Mev/c for the fol- 
lowing reason: The final-state interaction picture 
is valid only at low relative energies of the two 
pions where their attraction is large compared 

to other effects, such as details of the production 
mechanism and final-state 7-He’* interactions. 

In Fig. 4 we show the result of subtracting the 
phase-space volume from the data. The peak 
occurs at a mass value of about 300 Mev, and un- 
folding the resolution gives a width of about 25 Mev. 
This yields a lifetime of the same order as the 
interaction time, and the concept of a particle be- 
comes vague. It therefore seems doubtful that 
this is the vector meson of Nambu.*° We intend 
to settle this question independently of the reso- 
lution in a subsequent experiment. 

We cannot completely rule out the possibility 
that the anomaly is due to final-state interactions 
between the He® and one of the pions. The J= 3, 

J = 3, pion-nucleon state can occur in the H® case, 
but not as strongly as in the He* case. However, 
in the region of the anomaly we are below the en- 
ergy of the 3, 3 resonance. Also, in other ex- 
periments in which the final state consists of two 
pions and a nucleon, the energy spectrum of the 
nucleon is usually not strongly influenced by the 
8, 2 resonance." We can dismiss the S-wave 
pion-nucleus interaction because it is known to 
be small at these energies.” 

We have also considered the symmetrization 
of the wave function for the two pions. Details 
of the derivation are given in reference 3. We 
present here only the result, which is to multiply 
Eq. (5) by the factor B(R) =1+ exp{ -[(w/,)? - 4] 
x(R/2.15)?}, where R is the radius of interaction 
in units of l/u. The effect of the symmetrization 
on the shape of @, is small for all values of R 
and cannot reproduce the observed bump. We 
will therefore neglect B(R), although including 
it would decrease the value of the scattering 
length derived above. 

Another effect we mention only briefly is the 
He® wave function. Roughly speaking, the three 
nucleons stick together to form a He*® more easily 
when they have low relative energies because the 
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He® wave function has fewer high-momentum 
components. This effect favors low He* momenta, 
From both the H® and He® spectra, we conclude 
that this effect is small or else compensated for co 
by something else. a 

We plan to repeat the experiment at another 
laboratory-system angle. It should then be possi- 
ble to make the correct interpretation. 
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